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1.0  INTRODUCTION 


One  of  the  primary  funetions  of  the  Global  Oeean  Foreeast  System  (GOFS)  Version  3.0 
(V3.0)  (consisting  of  the  1/12°  global  HYbrid  Coordinate  Ocean  Model  (HYCOM)  that 
employs  the  Navy  Coupled  Ocean  Data  Assimilation  (NCODA),  Metzger  et  ah,  2008)  is  to 
provide  boundary  conditions  (BCs)  to  higher  horizontal  and  vertical  resolution  regional  nested 
ocean  models.  This  was  validated  in  the  GOFS  V3.0  Phase  II  Validation  Test  Report  (VTR) 
(Metzger  et  ah,  2010)  for  the  Relocatable  (Relo)  (or  Regional)  Navy  Coastal  Ocean  Model 
(NCOM)  (Rowley  et  ah,  2010)  configured  for  the  region  surrounding  Luzon  Strait  that  connects 
the  Pacific  Ocean  with  the  South  China  Sea.  The  Phase  II  VTR  determined  that  the  Luzon  Strait 
Relo  NCOM  hindcast  using  GOFS  V3.0  BCs  produced  superior  results  to  a  twin  hindcast  that 
used  BCs  from  GOFS  V2.6  (consisting  of  1/8°  global  NCOM/NCODA,  1/32°  Navy  Layered 
Ocean  Model  (NLOM)  and  1/8°  Modular  Ocean  Data  Assimilation  (MOD AS)). 

This  document  can  be  viewed  as  both  a  User’s  Manual  and  an  addendum  to  the  Phase  II 
boundary  condition  validation  work.  As  such,  the  dual  methodology  for  extracting  GOFS  V3.0 
BCs  is  documented  (Section  2).  Secondly,  it  describes  a  set  of  experiments  in  which  BCs 
extracted  from  the  real-time  GOFS  V3.0  have  been  remapped  to  40,  50  or  100  vertical  levels  and 
these  are  used  with  corresponding  versions  of  Relo  NCOM  configured  for  the  region  off  the  New 
Jersey  coast  (Section  3).  The  sensitivity  of  the  vertical  resolution  is  examined  with  regard  to 
error  analyses  of  temperature  (T)  vs.  depth  and  acoustical  proxy  measures  (i.e.  mixed  layer  depth 
(MLD),  sonic  layer  depth  (SLD),  below  layer  gradient  (BLG)  and  deep  sound  channel  (DSC)). 


Manuscript  approved  January  1 1,  2012. 


1 


2.0  METHODOLOGY  FOR  EXTRACTING  GOES  V3.0  BCs 


There  are  eurrently  two  sets  of  proeedures  for  extraeting  GOFS  V3.0  BCs  for  use  with 
Relo  NCOM.  Under  method  1,  as  part  of  the  real-time  system’s  runstream,  netCDF  files  are 
output  on  a  eonstant  .08°  latitude/longitude  grid  that  spans  the  globe  from  80°S-80°N.  (Currently 
these  netCDF  files  are  output  from  the  noweast  through  the  120-hr  foreeast  with  3-hourly 
temporal  frequeney.)  These  are  vertieally  interpolated  to  40  levels:  0,  2,  4,  6,  8,  10,  12,  15,  20, 
25,  30,  35,  40,  45,  50,  60,  70,  80,  90,  100,  125,  150,  200,  250,  300,  350,  400,  500,  600,  700,  800, 
900,  1000,  1250,  1500,  2000,  2500,  3000,  4000  and  5000  m.  A  sequenee  of  seripts  (defined  in 
Section  2.1)  is  then  executed  to  extract  the  BCs  for  a  user-defined  Relo  NCOM  domain.  Under 
method  2,  the  global  3D  HYbrid  Coordinate  Ocean  Model  (HYCOM)  archive  files  on  the  native 
vertical  hybrid  grid  (or  on  a  user-defined  subdomain)  can  be  extracted  and  saved  at  a  defined 
temporal  frequency.  A  different  sequence  of  scripts  (defined  in  Section  2.2)  is  then  executed  to 
extract  the  BCs  for  a  Relo  NCOM  domain.  An  advantage  of  the  second  method  is  that  the 
vertical  resolution  can  be  different  from  the  40  levels  defined  in  the  first  approach. 

The  scripts  defined  in  this  document  have  been  saved  using  the  Subversion  (SVN) 
configuration  management  software  and  can  be  found  at: 

setenv  SVNROOT  "https://www7320.nrlssc.navy.mil/svn/repos" 

${SVNROOT}/GOFSV3/trunk/postproc/BC_extraction/V1.0/script_cut 

${SVNROOT}/GOFSV3/trunk/postproc/BC_extraction/V1.0/script_glb 

$  { SVNROOT  }  /GOFS  V3/trunk/postproc/BC_extraction/V  1 . 0/ subregion 
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2.1  Method  1:  Procedure  to  extract  BCs  from  GOFS  V3.0  netCDF  files 


The  basic  idea  behind  both  methods  is  to  create  netCDF  files  with  the  same  naming 
convention  and  format  that  Relo  NCOM  uses  for  receiving  BCs  from  GOFS  V2.6.  Thus  it  is 
necessary  to  create  the  files  and  directory  structure  similar  to  the  hostnl  namelist  of  the  Relo 
NCOM  domain,  a  portion  of  which  is  noted  below: 


&hostnl 

host  ncpath  = 
host  odimens  = 
host  ohgrda  = 
host  ohgrdb  = 
host  ovgrdd  = 


' /scr/$ {user} /hycom/Nc_GLBuO . 08 
'/scr/${user} /hycom/GLBuO . 08 /input 0/odimens . D ' 
'/scr/${user} /hycom/GLBuO . 08/input0/ohgrd_l .A' , 
'/scr/${user} /hycom/GLBuO . 08/input0/ohgrd_l .B ' , 
'/scr/${user} /hycom/GLBuO . 08/input0/ovgrd_l . D '  , 


2.1.1  Script  set-up 


It  is  user’s  responsibility  to  create/copy/modify  the  files  and  to  define  the  environmental 
variables.  In  the  scripts  that  follow,  the  yellow  highlighted  areas  may  need  to  be  changed  by  the 
user. 


1.  The  file  /u/home/${user}/hycom_bin.env  contains  the  location  of  the  executables  and 
working  directories  required  for  extracting  GOFS  V3.0  BCs. 


/u/home/${user}/hycom_bin.env 

# 

#  -  This  file  defines  environmental  variables  that  are  used 

#  -  in  extracting  GOFS  V3 . 0  boundary  conditions  for  a  Relo 

#  -  NCOM  nested  domain 

# 

#  -  working  directory 

# 

setenv  WRK  dir  /scr/$ {user } /hycom 

# 

#  -  directory  for  scripts  to  extract  BCs: 

#  -  SCRPT  cut  is  the  methodology  for  using  HYCOM  archive  files 

#  -  that  are  a  subregion  of  the  entire  domain 

#  -  SCRPT  gib  is  the  methodology  for  using  the  constant  .08  deg 

#  -  lat/lon  netCDF  files 

# 

setenv  SCRPT_cut  /u/home/ $ { user } /hyc2ncom/script_cut 
setenv  SCRPT  gib  /u/home/ ${ user } /hyc2ncom/script  gib 
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# 

#  -  netCDF  operators  directory 

# 


setenv  NCO 
# 

#  -  executables 

#  _ 

setenv  AWBIN 

setenv  AWBIN2 
setenv  AWsub 
setenv  AWarc 
setenv  AWTOPO 
setenv  MODASbin 
setenv  PLTF 


/site/nco-3 . 9 . 8_64/bin/ 
and  utilities  directories 


/u /home /wallcraf/hycom/ ALL /bin| _ 

/u /home /wallcraf/hycom/ ALL /ncom/  src/| _ 

/u/home/wallcraf /hycom/ ALL/ subregion/ src| 
/u/home/wallcraf /hycom/ALL/archive/src| 
/u/home/wallcraf /hycom/ALL/topo  /  src| 
/u/home/  smedstad/bini 

navo 


# 

#  -  topography  directories:  note  GLBaO.08  is  the  native  Mercator- 

#  -  curvilinear  HYCOM  grid  whereas  GLBu0.08  is  the  constant  .08 

#  -  deg  lat/lon  grid  onto  which  the  netcdf  files  have  been 

#  -  interpolated 

# 

setenv  TOPOhy_dir  /u/home/ ${ user } /hycom/GLBaO . 08 /topo 
setenv  TOPOnc  dir  newton : /u/home/$ {user } /hycom/GLBuO . 08/topo 


2.  The  scripts  in  the  directory  /u/home/${user}/hyc2ncom/script_glb  are  used  on  the  GOFS  V3.0 
constant  .08°  lat/lon  netCDF  files,  i.e.  the  GLBuO.08  grid.  The  following  files  are  required; 


a.  /u/home/${user}/hyc2ncoin/script_glb/include.env  defines  more  environmental  variables, 
the  time  frame  over  which  to  extract  BCs,  the  GLBu0.08  array  size  and  the  GOFS  V3.0 
experiment  number,  in  this  example  908.  Typically,  dayl,  day2  and  maybe  taus  will  be  the  only 
variables  to  change  in  this  file. 


# ! /bin/csh 

# 

#  -  this  file  is  used  as  part  of  the  methodology  to  extract 

#  GOFS  V3 . 0  BCs  for  a  Relo  NCOM  domain.  In  this  case,  GOFS 

#  V3 . 0  input  files  are  in  netCDF  format  and  on  a  constant 

#  -  0.08  deg  lat/lon  grid 

# 

source  /u/home/$ {user } /hycom  bin.env 
setenv  SCRPT  $SCRPT_glb 

# 


setenv  AREA 
setenv  ncom  gr 
setenv  ARCHenv 
setenv  NC  area 
setenv  Topo_area 


GLBuO .08 

$WRK_dir/$AREA/inputO 
$WRK  dir/global  nc 
$ WRK_di r / Nc_$ ARE A 
$ARCHenv 
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# 

#  -  dayl  =  start  DTG,  day2  =  end  DTG,  taus  =  temporal  freq 

# 

Isetenv  dayl  201008^ 

Isetenv  day2  2010082TI 

Isetenv  taus  "000  003  006  009  012  015  018  021^ 

# 

#  -  array  dimensions  of  GOFS  V3 . 0  netcdf  files 

# 

setenv  NX  4500 
setenv  NY  2001 
# 

#  -  more  GOFS  V3 . 0  parameters: 

#  Tenv  =  topography  version  number 

#  REGenv  =  domain  name  for  the  constant  .08  deg  grid 

# 

Isetenv  Tenv  09| 
setenv  REGenv  GLBuO.08 
# 

#  -  expt  name  in  two  forms 

# 

setenv  Eenv  908 
setenv  Xenv  90 . 8 
# 

#  -  HRenv  =  analysis  time 

#  -  HR2env  =  nowcast  time 

# 

setenv  HRenv  18 
setenv  HR2env  00 

b.  /u/home/${user}/hyc2ncoin/script_glb/stepl.get_GOFSV3_ncdf.com  is  used  to  transfer  the 
GLBuO.08  netCDF  files  from  the  arehive  machine  (newton.navo.hpc.mil)  to  the  working 
directory. 


# ! /bin/csh  -x 

#PBS  -N  XXX 

#PBS  -j  oe 

#PBS  -o  XXX. log 

#PBS  -e  XXX. err 

#PBS  -1  walltime=6 : 00 : 00 

I# PBS  -A  NRLSS0181 

#PBS  -q  transfer 
#PBS  -1  select=l 
# 

set  echo 
set  time  =  1 
# 

cd 
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source  /u/home/$ {user } /hycom  bin.env 
source  $SCRPT  glb/include . env 
# 

#  -  location  of  GOFS  V3 . 0  netCDF  files  interpolated  to  constant 

#  -  0.08  deg  lat/lon  grid 

# 

Isetenv  V  /u/home/ooc/data/hycom/ $  { REGenv } /expt  $  {Xenv} /data/netcdf| 

setenv  S  $ARCHenv 

# 

if  (!  -e  ${S})  mkdir  -p  ${S} 
cd  ${S} 

|set  var  =  "ssh  ts3z  uv3^ 

Iset  tau  =  "000  003  006  009  012  015  018  021'^ 

# 

set  dstr  =  $dayl 
set  dend  =  $day2 
echo  $dstr  $dend 
# 

#  -  stage  files  on  newton 

# 

while  (  $dstr  <=  $dend) 
foreach  nm  ($var) 
foreach  tt  ($tau) 

rsh  newton  /opt/SUNWsamf s/bin/stage 
$ { V} /hycom_glb_$ {Eenv}_$ { dstr } 00_t$ { tt }_$ { nm} . nc 
end 
end 

set  dstr  =  ' $MODASbin/addndays  YYYYMMDD  $dstr  +1' 

end 

# 

set  dstr  =  $dayl 
set  dend  =  $day2 
# 

#  -  copy  files  from  newton 

# 

while  (  $dstr  <=  $dend  ) 
foreach  nm  ($var) 
foreach  tt  ($tau) 

/usr/bin/rcp 

newton :$ {V} /hycom_glb_$ { Eenv }_${ dstr } 00_t$ { tt }_$ {nm} . nc  .  & 
end 
end 
wait 

set  dstr  =  ' $MODASbin/addndays  YYYYMMDD  $dstr  +1' 

end 

c.  /u/home/${user}/hyc2ncom/script_glb/step2.manipulate_ncdf.com  is  used  to  manipulate 
the  GLBuO.08  netCDF  files  and  put  in  a  format  eonsistent  with  what  Relo  NCOM  expeets. 
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#  !  /bin/csh  -f 
#PBS  -N  XXX 
#PBS  -j  oe 
#PBS  -o  XXX. log 
#PBS  -e  XXX. err 

#PBS  -1  walltime=3 : 00 : 00 
I# PBS  -A  NRLSS0181 

#PBS  -q  standard 
#PBS  -1  select=l 

# 

set  echo 
set  time  =  1 
# 

#  -  this  script  manipulates  the  netcdf  files  and  modifies  to 

#  -  make  them  consistent  with  what  Relo  NCOM  wants 

# 

/bin/date 

cd 

source  /u/home/$ {user } /hycom  bin.env 
source  $SCRPT  glb/include . env 
# 

setenv  Nncks  $NCO/ncks 
setenv  WK  $NC  area 
mkdir  -p  $WK 

# 

setenv  ARCH  $ARCHenv 

# 

set  taus  =  '000  003  006  009  012  015  018  021' 
set  dstr  =  $dayl 
set  dend  =  $day2 
# 

set  var  =  "ssh  s3d  t3d  u3d  v3d" 
set  hnm  =  hycom  $Eenv 
# 

while ($dstr  <=  $dend) 
foreach  dd  ($taus) 

set  glfile  =  $ARCH/hycom  gib  $Eenv"  "$dstr"00  t"$dd 

/bin/cp  $glfile"  "ssh.nc  $WK/ssh  $hnm"  "$dstr"00  t"$dd"h".nc 
$Nncks  -0  -h  -v  water  temp  $glfile"  "ts3z.nc 
$WK/t3d  $hnm"  "$dstr"00  t"$dd"h".nc 

$Nncks  -0  -h  -v  salinity  $glfile"  "ts3z.nc 
$WK/s3d  $hnm"  "$dstr"00  t"$dd"h".nc 

$Nncks  -0  -h  -v  water  u  $glfile"  "uv3z.nc 
$WK/u3d  $hnm"  "$dstr"00  t"$dd"h".nc 

$Nncks  -0  -h  -v  water  v  $glfile"  "uv3z.nc 
$WK/v3d  $hnm"  "$dstr"00  t"$dd"h".nc 

foreach  name  ($var) 

set  file  =  $WK/$name"  "$hnm"  "$dstr"00  t"$dd"h".nc 
$NCO/ncwa  -0  -h  -a  time  $file  $file 
end 

echo  ncwa  $dd 

set  file  =  $hnm"  "$dstr"00  t"$dd"h".nc 
$Nncks  -0  -h  -v  surf  el  $WK/ssh  $file  $WK/ssh  $file 
$Nncks  -0  -h  -v  water  temp  $WK/t3d  $file  $WK/t3d  $file 
$Nncks  -0  -h  -v  salinity  $WK/s3d  $file  $WK/s3d  $file 
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$Nncks  -0  -h  -v  water  u  $WK/u3d  $file  $WK/u3d  $file 
$Nncks  -0  -h  -v  water  v  $WK/v3d  $file  $WK/v3d  $file 
end 

echo  $dstr  done 

set  dstr  =  ' $MODASbin/addndays  YYYYMMDD  $dstr  +1' 
end 
# 

/bin/date 

d.  /u/home/${user}/hyc2ncom/script_glb/get_ncgrid.com  is  used  to  transfer  from 

newton.navo.hpc.mil  the  files  associated  with  the  GOFS  V3.0  constant  .08  deg  lat/lon  grid. 

# ! /bin/csh 

#PBS  -N  XXX 

#PBS  -j  oe 

#PBS  -o  XXX. log 

#PBS  -e  XXX. err 

#PBS  -1  walltime=l : 00 : 00 

I# PBS  -A  NRLSS0181 

#PBS  -q  transfer 
#PBS  -1  select=l 
# 

#  -  script  to  transfer  from  newton  the  files  associated 

#  -  with  the  GOFS  V3 . 0  constant  .08  deg  lat/lon  grid 

# 

set  echo 
set  time  =  1 
# 

cd 

source  /u/home/$ {user } /hycom  bin.env 
source  $SCRPT  glb/include . env 
# 

setenv  V  /u/home/$ {user } /hycom/$AREA/inputO 
setenv  S  $ncom  gr 
# 

if  (!  -e  ${S})  mkdir  -p  ${S} 
cd  ${S} 

# 

#  -  stage  files  on  newton 

# 

set  file  =  "odimens.D  ohgrd  l.A  ohgrd  l.B  ovgrd  1 . D" 
foreach  ff  ($file) 

/usr/bin/rsh  newton  /opt/SUNWsamf s/bin/stage  ${V}/$ff 
end 
# 

#  -  copy  files  from  newton 

# 

foreach  ff  ($file) 

/usr/bin/rcp  newton :$ {V} /$ff  .  & 
end 
wait 
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2.1.2  Executing  the  scripts  to  extract  GOES  V3.0  BCs 


After  the  environmental  variables  are  set,  the  scripts  are  run  in  this  order. 

1.  stepl.get_GOFSV3_ncdf.com:  to  retrieve  the  GOFS  V3.0  netCDF  files  and  store  them  in 
SARCHenv 

2.  step2.manipulate_ncdf.com:  to  convert  the  GOFS  V3.0  netCDF  files  in  $ARCHenv  and 
store  the  new  format  in  $NC_area. 

3.  step3.get_GLBu0.08.com:  to  retrieve  the  $ncom_gr,  the  inputO  directory  of  hostnl 

4.  Specify  in  hostnl,  the  path  on  $ncom_gr  and  $  NC_area. 

5.  Run  Relo  NCOM. 

2.2  Method  2:  Procedure  to  extract  BCs  from  the  GOFS  V3.0  archive  files  on  the  native 
hybrid  vertical  grid 

Boundary  conditions  can  also  be  extracted  from  GOFS  V3.0  archive  files  on  the  native 
hybrid  vertical  grid.  This  procedure  is  more  cumbersome  but  has  the  advantage  of  allowing  the 
vertical  resolution  to  be  different  from  the  40  pre-defined  levels  used  in  the  GOFS  V3.0  netCDF 
files.  The  input  GOFS  V3.0  archive  files  can  cover  the  whole  domain,  or  just  a  subdomain 
(defined  in  the  archt.input  file  in  the  working  GOFS  V3.0  directory).  If  whole  domain  files  are 
used,  a  sequence  of  scripts  is  used  to  extract  a  portion  of  the  entire  domain  that  is  slightly  larger 
than  the  Relo  NCOM  domain.  These  smaller  files  are  then  easier  to  manage.  The  procedure  to 
process  the  GOFS  V3.0  subdomain  archive  files  is  detailed  in  the  appendix  and  this  must  be 
done  before  the  following  steps  are  performed.  Such  subdomain  files  would  have  to  be  created  as 
GOFS  V3.0  is  running. 

The  example  that  follows  is  set  up  for  a  Relo  NCOM  domain  that  covers  the  Luzon  Strait 
region.  The  longitude  and  latitude  ranges  are:  1 18. 825-126. 274°E  and  17.35-25. 153°N  and  it  is 
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configured  with  50  vertical  levels.  The  GOFS  V3.0  archive  files  must  be  slightly  larger  than  the 
Relo  NCOM  domain  and  this  example  will  use  1 18. 0-127. 0°E,  17.0-26.0°N.  The  user  must 
determine  the  number  of  GOFS  V3.0  gridpoints  in  the  west-east  (NX)  and  north-south  (NY) 
directions  and  this  is  done  as  follows: 


dl3n6  88>  /u/home/wallcraf /hycom/ALL/bin/hycom  lonlat2ij  118.0  17. 0| 
~/hycom/GLBaO  .  08/topo/reqional  .grid.  a| 

549  1721 _ 

dl3n6  89>  /u/home/wallcraf /hycom/ALL/bin/hycom  lonlat2iT  127.0  26. 0| 
-/hycom/GLBaO  .  08/topo/reqional  .grid.  a| 

661  1842 


So,  NX  =  661  -549+  1  =  113  and  NY  =  1842  -1721  +  1  =  122. 


1.  As  in  the  previous  methodology,  the  file  /u/home/${user}/hycom_bin.env  contains  the 
location  of  the  executables  and  working  directories  required  for  extracting  GOFS  V3.0  BCs. 

2.  The  scripts  in  the  directory  /u/home/${user}/hyc2ncoin/script_cut  are  used  on  the  GOFS 
V3.0  archive  files  on  the  native  hybrid  vertical  grid.  The  following  files  are  required: 

a.  /u/home/${user}/hyc2ncoin/script_cut/include.env  defines  more  environmental  variables, 
the  time  frame  over  which  to  extract  BCs,  the  array  size  and  the  experiment  number,  in  this 
example  908.  Typically  the  only  variables  that  change  in  this  file  are:  dayl/day2  -  the  start/end 
dates  and  NX/NY  -  the  longitude/latitude  dimensions  on  the  GLBa0.08  grid  (which  have  been 
determined  above). 


# ! /bin/csh 
# 

#  this  file  is  used  as  part  of  the  methodology  to  extract 

#  GOFS  V3 . 0  BCs  for  a  Relo  NCOM  domain.  In  this  case,  GOFS 

#  V3 . 0  input  files  are  a  subdomain  of  the  full  model  grid 

#  and  on  the  hybrid  vertical  grid,  i.e.  not  z-levels 

# 

source  /u/home/$ {user } /hycom  bin.env 
setenv  SCRPT  $SCRPT  cut  “ 
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# 

#  -  array  dimensions  of  the  GOFS  V3 . 0  subdomain 

# 


setenv 

NX 

113 

setenv 

NY 

122 

# 

#  -  environmental  variable  to  define  how  the  Relo  NCOM  vertical 

#  -  is  configured.  When  set  to  1,  use  pre-defined  depths,  when  set 

#  -  to  0  use  the  Relo  NCOM  algorithm 

# 

setenv  MAKE_VGRID  1 

# 

#  -  dayl  =  start  DTG,  day2  =  end  DTG 

#  _ 

setenv  dayl  20100702 
setenv  day2  20100706 

setenv  year  'echo  $dayl  |  awk  '{print  substr ($1, 1, 4) } ' ' 

# 

#  -  directories 

# 

setenv  TMP  $WRK  dir/tmp 

setenv  AREA  LuzonStrait  cut 

setenv  ncom  gr  $WRK  dir/$AREA/inputO 
setenv  Topo  area  /u/home/$ {user } /hyc2ncom/Topo  $AREA 
setenv  ARCHenv  $WRK  dir/global 
setenv  TOPO  dir  $TOPOhy  dir 
setenv  ARCH  area  $WRK  dir/Arch  $AREA 
setenv  NC  area  $WRK  dir/Nc  $AREA 
#  “ 

if  ( !  -e  $Topo  area)  mkdir  -p  $Topo  area 

if  ( !  -e  $ARCH  area)  mkdir  -p  $ARCH  area 

if  ( !  -e  $NC  area)  mkdir  -p  $NC  area 
#  “ 

#  more  GOFS  V3 . 0  parameters: 

#  -  Tenv  =  topography  version  number 

#  -  REGenv  =  domain  name  for  GOFS  V3 . 0 

# 

setenv  Tenv  09 
setenv  REGenv  GLBaO.08 
# 

#  -  expt  number 

# 

setenv  Eenv  908 

# 

#  -  HRenv  =  analysis  time 

#  -  HR2env  =  nowcast  time 

# 

setenv  HRenv  18 
setenv  HR2env  00 
# 

#  -  define  Julian  dates 

# 

if  (  $PLTF  ==  navo)  then 

set  mm  =  'echo  $dayl  |  awk  '{print  substr  ( $1 , 5 , 2 )  } ' ' 

set  dd  =  'echo  $dayl  |  awk  '{print  substr  ( $1 , 7 , 2 )  } ' ' 

set  jj  =  ' csh  $MODASbin/ymd2doy . csh  $year  $mm  $dd' 
setenv  jull  'echo  $jj  |  awk  '{print  substr ( $1 , 6 , 3 ) } ' ' 


11 


set  mm  =  'echo  $day2  |  awk  '{print  substr ( $1 , 5 , 2 ) } ' ' 
set  dd  =  'echo  $day2  |  awk  '{print  substr ( $1 , 7 , 2 ) } ' ' 
csh  $M0DASbin/ymd2doy . csh  $year  $mm  $dd 
set  jj  =  'csh  $M0DASbin/ymd2doy . csh  $year  $mm  $dd' 
setenv  jul2  'echo  $jj  |  awk  '{print  substr ( $1 , 6 , 3 ) } ' ' 
else 

setenv  jull  ' $M0DASbin/idtg2 jul  $dayl ' 
setenv  jul2  ' $M0DASbin/idtg2 jul  $day2 ' 
endif 

b.  /u/home/${user}/hyc2ncoin/script_cut/stepl.get_GOFSV3_arch.com  will  transfer  from 
newton.navo.hpc.mil  the  tarballs  that  contains  the  best  GOFS  V3.0  analysis  for  each  day  of  a 
specific  month  and  year. 


# ! /bin/csh 

# 

#PBS  -N  XXX 

#PBS  -j  oe 

#PBS  -o  XXX. log 

#PBS  -e  XXX. err 

#PBS  -1  walltime=6 : 00 : 00 

I# PBS  -A  NRLSS0181 

#PBS  -q  transfer 
#PBS  -1  select=l 
# 

set  echo 
set  time  =  1 
# 

#  -  this  script  will  transfer  from  newton  the  tarballs  that  contains 

#  -  the  best  GOFS  V3 . 0  analysis  for  each  day  of  a  specific  month  and 

#  - year 

# 

cd 

source  ~/hycom  bin.env 
source  $SCRPT  cut/include . env 
# 

setenv  R  GLBaO.08 
setenv  E  908 
setenv  X  90.8 
setenv  Y  2010 

setenv  M  07 _ 

Isetenv  V  /u/home/ooc/data/hycom/ $  { R} /expt  ${X}/data| 

setenv  S  / scr/ $ { user } /hycom/global 

# 

if  ( !  -e  $ { S  } )  then 
mkdir  -p  ${S} 
endif 
cd  ${S} 

# 

#  -  list  all  best  analysis  tarballs  for  a  given  month  into  a  file; 

#  -  sort  through  and  find  the  best  analysis 

# 

#  -  the  first  day  of  the  month  is  special 

# 
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rsh  newton  "cd  ${V};  Is  ${E}  archv  ${Y}??????  $ { Y} $ {M} 0100*gz"  >!  11 
rsh  newton  "cd  ${V};  Is  ${E}  archv  ${Y}??????  ${Y}${M}*gz"  >!  12 

cat  11  12  >!  listl 

# 

set  L  =  'wc  -1  listl  I  awk  ' {printf ( "%d" ,  $1-1)}'' 
tail  -n  ${L}  listl  >!  Iist2 
# 

cat  list2  I  sed  -e  's/_/  /g'  -e  's/\./  /g'  |  sort  -nr  -k  4  |  uniq  -f  3 
I  sort  -n  -k  4  |  sed  -e  's/  gz/.gz/g'  -e  ' s/  tar/.tar/g'  -e  's/  /_/g' 

>!  listl 

# 

set  L  =  'wc  -1  listl  I  awk  ' {printf ( "%d" ,  $1)}'' 

# 

cat  listl 

# 

#  -  copy  files  from  newton 

# 

@  N  =  1 

while  (${N}  <=  ${L}) 

set  F  =  'head  -n  ${N}  listl  |  tail  -n  1' 
rep  newton :$ {V} /${ F}  .  & 

@  N  ++ 
end 
wait 
# 

#  -  untar/uncompress  tarballs 

# 

@  N  =  1 

while  (${N}  <=  ${L}) 

set  F  =  'head  -n  ${N}  listl  |  tail  -n  1' 

/site/unsupported/bin/gtar  --f ormat=posix  -xvzf  ${F}  & 

@  N  ++ 
end 
wait 
# 

#  -  delete  tarballs 

# 

@  N  =  1 

while  (${N}  <=  ${L}) 

set  F  =  'head  -n  ${N}  listl  |  tail  -n  1' 

/bin/rm  -f  ${F} 

@  N  ++ 
end 
# 

/bin/rm  -f  list[12]  1[12] 

c.  /u/home/${user}/hyc2ncom/script_cut/step2.get_GLBa0.08.com  will  link  or  transfer  from 
newton.navo.lipc.mil  the  needed  topography  and  latitude/longitude  definition  files. 


# ! /bin/esh 

# 

#PBS  -N  XXX 
#PBS  -j  oe 
#PBS  -o  XXX. log 
#PBS  -e  XXX. err 
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#PBS  -1  walltime=6 : 00 : 00 
I# PBS  -A  NRLSS0181 

#PBS  -q  transfer 
#PBS  -1  select=l 
# 

set  echo 
set  time  =  1 
# 

cd 

source  /u/home/$ {user } /hycom  bin.env 
source  $SCRPT  cut/include . env 

setenv  V  /u/home/$ {user} /hycom/$ {REGenv} /topo 
setenv  S  $ARCHenv 
setenv  LINK  1 
# 

if  (!  -e  ${S})  mkdir  -p  ${S} 
cd  ${S} 

# 

if  ($LINK  ==  1)  then 
/bin/rm  regional.* 

In  -sf  $ {V} /regional . grid . a  . 

In  -sf  $ {V} /regional . grid . b  . 

In  -sf  ${V}/depth  ${REGenv}  ${Tenv}.b  regional . depth . b 
In  -sf  ${V}/depth  ${REGenv}  ${Tenv}.a  regional . depth . a 
else 

rsh  newton  /opt/SUNWsamf s/bin/stage  $ {V} /regional . grid. [ab] 
rep  newton :$ {V} /regional . grid . a  .  & 
rep  newton :$ {V} /regional . grid . b  .  & 

rsh  newton  /opt/SUNWsamf s/bin/stage  ${V}/depth  ${REGenv}  ${Tenv}* 
rep  newton :$ {V} /depth  ${REGenv}  ${Tenv}.a  regional . depth . a 
rep  newton :$ {V} /depth  ${REGenv}  ${Tenv}.b  regional . depth . b 
endif 
wait 


3.  Next  the  Relo  NCOM  horizontal  grid  must  be  eonfigured.  It  is  assumed  the  user  has  suffieient 
knowledge  to  do  this.  Edit  /u/home/${user}/hyc2ncom/script_cut/relo.nl.  Here,  variables  m  = 
$NX  and  n  =  $NY,  rlon/rlat  are  the  starting  longitude/latitude,  dmax  is  the  deepest  level  and 
the  first  element  of  lo  defines  the  number  of  levels.  NOTE:  STEPS  3-7  ONLY  HAVE  TO  BE 
DONE  ONE  TIME. 


Sdsetnl 

/ 

Sgridnl 

nnest  = 

1, 

nproj  = 

5, 

m  = 

113 

f 

!  longitude 

dimensions 

n  = 

122 

f 

!  latitude 

dimensions 

rlon  = 

118 

!  starting 

longitude 

rlat  = 

17  . 

!  starting 

latitude 

m  must  be  before  n 
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iref  =  1, 
jref  =  1, 
ii  =1, 

j  j  =1^ 

delx  =  0.0833333, 
dely  =  0.0833333, 

/ 

&hostnl 

/ 

&oanl 

/ 

Somnlof f 

/ 

Sparmlst 

/ 

&rlx3nl 


/ 

&setupl 


bathyfile  = 

' /u/home/ooc/models/relo/relo 

etc/dbdb2  v30. 

dat '  ,1 

dmax  = 

1-55001.  , 

dmin  = 

-10., 

dztop  = 

.5, 

qdem  dir  = 

' /u/home/ rowley/usr/relo/etc//ncoda/gdem3s ' , 

Qdemfile  = 

' /u/home/rowley/usr 

/relo/etc/cidem3  tsO.dat', 

initialtide  = 

. false . , 

lo  = 

Mr 

2, 

50, 

50, 

50, 

50, 

50, 

Iso  = 

2, 

35, 

35, 

35, 

35, 

35, 

35, 

nobmaxo  = 

4000, 

4000, 

4000, 

1000, 

1000, 

1000, 

1000, 

nqo  = 

2, 

2, 

2, 

2, 

2, 

2, 

2, 

nrivo  = 

200, 

50, 

50, 

50, 

50, 

50, 

50, 

nro  = 

2, 

2, 

2, 

2, 

2, 

2, 

2, 

ntco  = 

8, 

8, 

8, 

8, 

8, 

8, 

8, 

ntypo  = 

Ir 

Ir 

Ir 

1, 

Ir 

Ir 

Ir 

riverfile  = 

' /net/dynamic/export/data/rowley/models/relo/relo  1 . l/etc/rivers6 . dat ' 


startatrest  =  .false., 

_ tidefile  = _ 

' /net/dynamic/export/data/rowley/models/relo/relo_l . l/etc/tide_eqb . dat| 


writeinit 

writeosstf 

writeotsf 


. true . , 

. false . , 
. false . , 


/ 


&sf Ixnl 

/ 
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a.  There  are  two  options  for  creating  the  Relo  NCOM  vertical  grid  (ovgrid)  and  these  are 
controlled  by  the  environmental  variable  MAKE  VGRID  found  in  include.env.  The  scripts 
/u/home/${user}/hyc2ncoin/script_cut/step3.make_reloncom.*  define  the  procedure.  If 
MAKE  VGRID  =  I,  (i.e.  pre-defined  depths),  an  ASCII  file  (zin.dat)  with  the  appropriate 
number  of  depths  and  levels  must  exist  in  $SCRPT_cut/bin. 


/u/home/${user}/hyc2ncoin/script_cut/bin/zin.dat 


50 

*  number 

of 

levels 

0.000 

*  z 

*  = 

sample 

depth 

0.500 

*  z 

*  = 

sample 

depth 

1.083 

*  z 

*  = 

sample 

depth 

1.762 

*  z 

*  = 

sample 

depth 

2.554 

*  z 

*  = 

sample 

depth 

3 .477 

*  z 

*  = 

sample 

depth 

4.552 

*  z 

*  = 

sample 

depth 

5.806 

*  z 

*  = 

sample 

depth 

7.268 

*  z 

*  = 

sample 

depth 

8 . 971 

*  z 

*  = 

sample 

depth 

10 . 957 

*  z 

*  = 

sample 

depth 

13.271 

*  z 

*  = 

sample 

depth 

15.968 

*  z 

*  = 

sample 

depth 

19.112 

*  z 

*  = 

sample 

depth 

22 .777 

*  z 

*  = 

sample 

depth 

27 .049 

*  z 

*  = 

sample 

depth 

32 . 027 

*  z 

*  = 

sample 

depth 

37 . 831 

*  z 

*  = 

sample 

depth 

44.595 

*  z 

*  = 

sample 

depth 

52 .479 

*  z 

*  = 

sample 

depth 

61.669 

*  z 

*  = 

sample 

depth 

72.380 

*  z 

*  = 

sample 

depth 

84 .865 

*  z 

*  = 

sample 

depth 

99.418 

*  z 

*  = 

sample 

depth 

116.380 

*  z 

*  = 

sample 

depth 

136.151 

*  z 

*  = 

sample 

depth 

159.195 

*  z 

*  = 

sample 

depth 

186.055 

*  z 

*  = 

sample 

depth 

217.363 

*  z 

*  = 

sample 

depth 

253 . 855 

*  z 

*  = 

sample 

depth 

296.390 

*  z 

*  = 

sample 

depth 

345.968 

*  z 

*  = 

sample 

depth 

403 .755 

*  z 

*  = 

sample 

depth 

471 .110 

*  z 

*  = 

sample 

depth 

549.619 

*  z 

*  = 

sample 

depth 

641 .128 

*  z 

*  = 

sample 

depth 

747 .789 

*  z 

*  = 

sample 

depth 

872 .111 

*  z 

*  = 

sample 

depth 

1017 .019 

*  z 

*  = 

sample 

depth 

1185 . 922 

*  z 

*  = 

sample 

depth 

1382 .793 

*  z 

*  = 

sample 

depth 

16 


1612.263 
1879.729 
2191 .483 
2554 .860 
2978.406 
3472 .000 
4047 .000 
4718.000 
5500 .000 


'  z 
'  z 
'  z 
'  z 
'  z 
'  z 
'  z 
'  z 


z 


=  sample  depth 
=  sample  depth 
=  sample  depth 
=  sample  depth 
=  sample  depth 
=  sample  depth 
=  sample  depth 
=  sample  depth 
=  sample  depth 


If  MAKE  VGRID  =  0,  the  Relo  NCOM  algorithm  will  be  used  to  define  the  depths  and  levels. 

/u/home/${user}/hyc2ncoin/script_cut/step3.make_reloncom.com 

#!/bin/esh 

source  /u/home/$  {user}/hycom_bin.env 
source  $SCRPT_cut/include.env 

$SCRPT/step3.make_reloncom.s 

/u/home/${user}/hyc2nconi/script_cut/step3.make_reloncom.s 

# ! /bin/sh 


[export  RELO=/u/home/rowley/usr/relo| 

export  BINDIR=$RELO/bin 
export  JOBDIR=$RELO/ script 
export  ETCDIR=$RELO/etc 
export  DEFDIR=$RELO/def ault 
export  STATIC=$WRK_dir 
export  REGION=$AREA 
mkdir  -p  $STATIC/$REGION 
mkdir  -p  $STATIC/$REGION/inputO 

cd  $STATIC/$REGION 

export  NCOM  OHGRD  lB=input0/ohgrd  l.B 
export  NCOM  OHGRD  lA=input0/ohgrd  l.A 
export  NCOM  DIMEN  0D=input0/odimens . D 
export  NCOM  OVGRD  lD=input0/ovgrd  l.D 
export  NCOM  OZOUT  lD=input0/ozout  l.D 

#ovgrid  from  relonl 

if  [  $MAKE_VGRID  -eg  0  ] ;  then 

$BINDIR/ncom  config.xc  -dovgrd  relo.nl  >input0/ovgrid. log 
else 

csh  $SCRPT/make  ovgrid.com 
fi 

#  make  the  ohgrid  for  Relo  NCOM 
/bin/cp  $SCRPT/relo.nl  $STATIC/$REGION 
$BINDIR/ncom  config.xc  -dohgrd  relo.nl 

cd  $SCRPT 

/u/home/${user}/hyc2ncoin/script_cut/make_ovgrid.com 
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# ! /bin/csh 

# 

source  /u/home/$ {user } /hycom  bin.env 
source  $SCRPT  cut/include . env 
# 

cd  $SCRPT/bin 

# 

/bin/rm  odimens.D 
/bin/rm  ovgrd  l.D 

In  -s  $ncom  gr/odimens.D  odimens.D 
In  -s  $ncom  gr/ovgrd  l.D  ovgrd  l.D 
# 

$SCRPT/bin/z2ovgrid. exe  <<E-o-D 
Sinputs 
Nx=$NX 
Ny=$NY 

/ 

E-o-D 
cd  $SCRPT 


b.  /u/home/${user}/hyc2ncoin/script_cut/step4.make_regiongrid.com  creates  a  HYCOM 
version  of  the  Relo  NCOM  grid  based  on  the  NCOM  longitude  and  latitude  fields. 


# ! /bin/ csh 

# 

set  echo 

# 

cd 

source  /u/home/$ {user } /hycom  bin.env 
source  $SCRPT  cut/include . env 
# 

#  -  a  HYCOM  version  of  a  NCOM  grid. 

#  -  based  on  the  NCOM  longitude  and  latitude  fields. 

# 

setenv  N  $ncom  gr 
setenv  IDM  $NX  “ 
setenv  JDM  $NY 
# 

cd  $Topo_area 

/bin/rm  -f  regional . lonlat . a+  regional . lonlat . a 

# 

$AWBIN/raw2hycom  $N/ohgrd_l.A  ${IDM}  ${JDM} 

regional . lonlat . a+ 

$AWBIN/hycom  extract  regional . lonlat . a+  ${IDM}  ${JDM}  2111 

regional . lonlat . a 

/bin/rm  -f  regional . lonlat . a+ 

# 

#  -  temporary  regional . grid. b 

# 

cat  >!  regional . grid . b  <<E-o-D 

${IDM}  ' idm  '  =  longitudinal  array  size 

${JDM}  'jdm  '  =  latitudinal  array  size 

E-o-D 

setenv  FOR051A  regional . lonlat . a 
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setenv  FOR061A  fort.61A 

$AWTOPO/grid  lonlat  2d  <<E-o-D 

${IDM}  ' idm  '  =  longitudinal  array  size 

${JDM}  'jdm  '  =  latitudinal  array  size 

E-o-D 

mv  fort. 61  regional . grid. b 
mv  fort.GlA  regional . grid. a 

cd  $SCRPT 

c.  /u/home/${user}/hyc2ncoin/script_cut/step5.make_gmapi.com  prepares  the  weight  function 

of  the  interpolation  from  2D  native  grid  to  2D  regular  grid. 

#!bin/csh  -f 
set  echo 
set  time=l 
# 

cd 

source  /u/home/$ {user } /hycom  bin.env 
source  $SCRPT  cut/include . env 
# 

#  -  form  subregion  grid  array  index  map  file,  GLBaO.08  to  the 

#  -  Luzon  Strait  cutout  domain 

# 

setenv  R  $Topo_area 
file  $ {R} /regional . grid. a 
# 

mkdir  -p  $TMP 
cd  $TMP 

if (  ! (-e  regional . grid. a) )  then 

/bin/cp  $TOPO  dir/regional . grid. a  . 

/bin/cp  $TOPO  dir/regional . grid. b  . 
endif 

if (  ! (-e  regional . depth . a) )  then 

/bin/cp  $TOPO  dir/regional . depth . a  . 

/bin/cp  $TOPO  dir/regional . depth .b  . 
endif 

touch  ${ R} /regional . gmapi  $REGenv.a 
/bin/rm  ${ R} /regional . gmapi  $REGenv. [ab] 

# 

$AWsub/isuba  gmapi  <<E-o-D 

${R}/regional.grid.a 

${ R} /regional . gmapi  $REGenv.a 

$REGenv  (4500x3298)  to  $AREA  (maxinc=25) 

$NX  'idm  '  =  longitudinal  array  size  of  subregion 

$NY  'jdm  '  =  latitudinal  array  size 

25  'maxinc'  =  maximum  input  array  index  jump  on  target  grid 

E-o-D 

# 

cd  $WRK  dir 

d.  /u/home/${user}/hyc2ncom/script_cut/step6.depth_GLBu0.08.com  forms  a  subregion 
bathymetry  file,  GLBaO.08  to  GLBuO.08. 
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# ! /bin/csh 

# 

set  echo 

# 

source  /u/home/$ {user } /hycom  bin.env 
source  $SCRPT  cut/include . env 
# 

#  -  form  subregion  bathymetry  file,  GLBaO.08  to  Luzon  Strait  cutout 

domain 

# 

setenv  R  $Topo_area 

# 

mkdir  -p  $TMP 
cd  $TMP 
# 

if (  ! (-e  regional . grid. a) )  then 

/bin/cp  $TOPO  dir/regional . grid. a  . 

/bin/cp  $TOPO  dir/regional . grid. b  . 
endif 

if (  ! (-e  regional . depth . a) )  then 

/bin/cp  $TOPO  dir/regional . depth . a  . 

/bin/cp  $TOPO  dir/regional . depth .b  . 
endif 
# 

touch  $ { R} /depth_$ { AREA}_$ { Tenv} . [ab] 

/bin/rm  $ {R} /depth_$ {AREA}_$ {Tenv} . [ab] 

# 

$AWsub/isuba  topog  <<E-o-D 
${ R} /regional . gmapi  GLBaO.OS.a 
$TOPO  dir/depth  ${REGenv}  ${Tenv}.b 
$ { R} / depth_$ { AREA}_$ { Tenv} . b 

depth  GLBaO.08  subregioned  to  GLBu0.08  via  isuba  topog 
$NX  'idm  '  =  longitudinal  array  size  of  subregion 

$NY  'jdm  '  =  latitudinal  array  size 

E-o-D 

e.  /u/home/${user}/hyc2ncoin/script_cut/step7.reg_2_ncomgr.com  remakes  the  ohgrd  Relo 
NCOM  file  with  the  GOFS  V3.0  bathymetry. 


set  echo 

# 

#  -  form  NCOM  ohgrd  file  for  GLBaO.08 

# 

source  /u/home/$ {user } /hycom  bin.env 
source  $SCRPT  cut/include . env 
setenv  D  $Topo_area 
setenv  R  $ncom  gr 
# 

/bin/rm  regional . grid . * 

In  -sf  $D/regional . grid. a  . 

In  -sf  $D/regional . grid. b  . 

# 

/bin/rm  regional . depth . * 

In  -s  $D/depth  $AREA"  "$Tenv.a  regional . depth . a 
In  -s  $D/depth  $AREA"  "$Tenv.b  regional . depth . b 
# 
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setenv  NCOM_DIMEN_OD  $R/odimens.D 
setenv  NC0M_0VGRD_1D  $R/ovgrd_l.D 
# 

setenv  NC0M_0HGRD_1A  $R/ohgrd_l.A 
setenv  NC0M_0HGRD_1B  $R/ohgrd_l.B 
# 

touch  $NC0M_0HGRD_1A  $NC0M_0HGRD_1B 
touch  regional . dncom . b 
/bin/rm  regional . dncom. [ab] 

/bin/mv  $R/ohgrd  l.A  $R/ohgrd  l.A.relo 
/bin/mv  $R/ohgrd  l.B  $R/ohgrd  l.B.relo 
# 

$AWBIN2/grid2ncom  <<E-o-D 

1  'ilst  '  =  1st  hycom  i-point  on  ncom  grid 

1  'jlst  '  =  1st  hycom  j-point  on  ncom  grid 

E-o-D 

# 

/bin/rm  regional . grid. [ab]  regional . depth . [ab] 


4.  Convert  the  GOFS  V3.0  archive  files  to  Relo  NCOM  netCDF  files.  These  scripts  are  computer 
intensive. 


a.  /u/home/${user}/hyc2ncoin/script_cut/step8.arch_2_ncdf.com  forms  interpolated  subregion 
archive  files,  GLBaO.08  to  GLBuO.08. 


# ! /bin/csh 

#PBS  -N  XXX 

#PBS  -j  oe 

#PBS  -o  XXX. log 

#PBS  -e  XXX. log 

#PBS  -1  walltime=2 : 00 : 00 

I# PBS  -A  NRLSS0181 

#PBS  -q  standard 

#PBS  -1  select=l : ncpus=l 

# 

# 

set  echo 
set  time=l 
# 

set  OS='uname' 
switch  ($OS) 
case  ' Linux ' : 

which  aprun 

if  (!  $status)  then 

set  APRUN= ' aprun  -n  1  -m  15g  ' _ 

|set  SRC= /u /home /wallcraf /hycom/ ALLcnll 

else 

set  APRUN='' _ 

|set  SRC= /u /home /wallcraf /hycom/ ALL| 

endif 
breaksw 
case  'AIX' 
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set  APRUN=' ' 

|set  SRC=/u/home/wallcraf /hycom/ALiJ 

breaksw 
default : 

set  APRUN='  ' _ 

|set  SRC=/u/home/wallcraf /hycom/ALlJ 

endsw 

# 

date 

# 


cd 

source  /u/home/$ {user } /hycom  bin.env 
source  $SCRPT  cut/include . env 
# 

#  -  form  interpolated  subregion  archive  files,  GLBaO.08  to 

#  -  Luzon  Strait  cutout  domain 

# 

#  -  R  is  the  original  region 

#  -  U  is  the  target  region 

#  -  D  is  the  location  of  the  original  archive  files. 

#  -  N  is  the  location  of  the  subregion  archive  files. 

#  -  E,y,d,h  select  the  archive  files. 

#  -  T  is  topography  number. 

# 

setenv  T  $Tenv 
set  REG=$REGenv 
setenv  U  $AREA 
set  E=$Eenv 

setenv  X  'echo  ${E}  |  awk  ' {printf ( " %04 . If " ,  $1*0.1)}'' 

set  EXPT=expt_${X} 

# 

setenv  idtg  $dayl 
setenv  idtgtod  $day2 

#  -  analysis  time 

setenv  HR  $HRenv 

#  -  nowcast  time 

setenv  HR2  $HR2env 

# 

setenv  INP  $ARCHenv 
setenv  OUT  $ARCH  area 
mkdir  -p  “${OUT} 

Ifs  setstripe  -d  ${OUT} 

Ifs  setstripe  ${OUT}  1048576  -1  8 

# 


cd  ${INP} 

# 

if (  ! (-e  ${ INP } /regional . grid. a) )  then 

/bin/cp  $TOPO  dir/regional .  grid .  a  ${ INP } /regional . grid . a 

/bin/cp  $TOPO  dir/regional . grid . b  ${ INP } /regional . grid . b 

endif 

if  (  !  (-e  ${ INP } /regional . depth . a) )  then 

/bin/cp  $TOPO_dir/depth_$ {REG}_$ {T } . a  ${ INP } /regional . depth . a 

/bin/cp  $TOPO_dir/depth_$ { REG}_$ { T } . b  ${ INP } /regional . depth . b 

endif 

# 

if (  ! (-e  $ {OUT } /regional . gmap . a) )  then 

/bin/cp  $Topo_area/ regional . gmapi_$ {REG } . a  $ { OUT } / regional . gmap . a 
/bin/cp  $Topo  area/regional . gmapi  ${REG}.b  ${ OUT } /regional . gmap . b 
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endif 

if (  ! (-e  $ {OUT } /regional . grid. a) )  then 
/bin/cp  $Topo_area/ regional . grid,  a 
/bin/cp  $Topo  area/regional . grid. b 
endif 

if (  ! (-e  $ {OUT } /regional . depth . a) )  then 
/bin/cp  $Topo_area/depth_$ {U}_$ {T } .a 
/bin/cp  $Topo_area/depth_$ {U}_$ {T } . b 
endif 
# 

set  jul=$jull 
@  jul  =  $jul  -  1 
while  ($jul  <=  $jul2) 

# 


$ { OUT } / regional . grid . a 
$ { OUT } / regional . grid . b 


$ { OUT } / regional . depth . a 
$ { OUT } / regional . depth . b 


if  (-e  $ {OUT } /archv_$AREA. $ { year }_$ { jul }_$ {HR2 } . a)  then 
/bin/rm  ${OUT}/archv  $AREA . $ { year }  ${jul}  $ { HR2 } . * 
endif 


${APRUN}  $ { AWsub } /isubaregion  <<E-o-D 
$ { OUT } /regional . grid . a 
$ { OUT } / regional . gmap . a 
$ { OUT } /regional . depth . a 
${INP}/regional. depth .  a 

$ { INP } /$ {E }_archv . $ { year }_$ { jul }_$ {HR2 } . a 
$ { OUT } / archv_$AREA . $ { year }_$ { j  ul }_$ { HR2 } . a 
${REG}  interpolated  to  ${U} 


$NX 

'idm 

=  target 

longitudinal  array 

size 

$NY 

'  jdm 

=  target 

latitudinal  array 

size 

1 

model ) 

'  icef Ig ' 

=  ice  in 

output  archive  flag 

(0=none, 1 

0 

E-o-D 

@  jul  = 
end 
# 

date 

' smooth ' 

$jul  +  1 

=  smooth 

interface  depths 

o 

II 

VI 

h- ' 

II 

v 

b.  /u/home/${user}/hyc2ncoin/script_cut/step9.ncdf_vertinterp.com.*  remaps  the  arehive 
files  to  the  defined  vertical  levels  and  creates  netCDF  files  for  Relo  NCOM. 


step9.ncdf_vertinterp.com 

# ! /bin/csh 

#PBS  -N  XXX 

#PBS  -j  oe 

#PBS  -e  XXX. err 

#PBS  -o  XXX. log 

#PBS  -1  walltime=0 : 2 0 : 00 

#PBS  -A  NRLSS018 

#PBS  -q  standard 

#PBS  -1  select=l : ncpus=l 

# 

set  echo 
set  time  =  1 
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# 

cd 

source  /u/home/$ {user } /hycom  bin.env 
source  $SCRPT  cut/include . env 
# 

$SCRPT/step9 . ncdf  vertinterp.s 

step9.ncdf_vertinterp.s 


# ! /bin/sh 

# 

set  echo 
time  =  1 
# 

#  -  Script  to  "cutout"  a  hycom  archive  file  to  a  netcdf  file 

#  -  first  do  the  2d  files  and  then  the  3D  files 

#  -  3D  files  are  extracted  from  hybrid  into  Z-coordinates 

# 

ncdf 2d  0  { 

# 

#  -  extract  2D  fields  from  a  single  HYCOM  archive  file 

# 

CDF023=${NCDIR}/ssh"_"${DOMAIN} "_"${E} "_"${dtg}${tau} . nc 
mkdir  -p  $NCDIR 
# 

touch  $CDF023 
/bin/rm  $CDF023 
# 

export  CDF  TITLE="HYCOM  $DOMAIN  $EXPT" _ 

[export  CDF  INST="Naval  Research  Laboratory"! 

export  CDF023=$CDF023 
$AWarc/archv2ncdf 2d  <<  EOF 
$ { HYCOMINPUTFILE } 

MERSEA 


000 

'iexpt  ' 

= 

experiment  number  xlO  (000=from  archive 

file) 

3 

' yrf lag ' 

= 

days  in  year  flag  (0=360J16, 1=366J16, 2= 

366J01, 3 

actual ) 

$NX 

'  idm  ' 

= 

longitudinal  array  size 

$NY 

'  jdm 

= 

latitudinal  array  size 

1 

'  kdm  ' 

= 

number  of  layers 

34.0 

'  thbase ' 

= 

reference  density  (sigma  units) 

0 

' smooth ' 

= 

smooth  fields  before  plotting  (0=F,1=T) 

0 

'mthin  ' 

= 

mask  thin  layers  from  plots  (0=F,1=T) 

1 

'  iorign ' 

= 

i-origin  of  plotted  subregion 

1 

'  j  orign ' 

= 

j -origin  of  plotted  subregion 

0 

' i dmp  ' 

= 

i-extent  of  plotted  subregion  (<=idm;  0 

implies 

idm) 

0 

' j  dmp  ' 

= 

j -extent  of  plotted  subregion  (<=jdm;  0 

implies 

j  dm) 

0 

'botio  ' 

= 

bathymetry  I/O  unit  (0  no  I/O) 

0 

'flxio  ' 

= 

surf,  heat  flux  I/O  unit  (0  no  I/O) 

0 

'empio  ' 

= 

surf,  evap-pcip  I/O  unit  (0  no  I/O) 

0 

'ttrio  ' 

= 

surf,  temp  trend  I/O  unit  (0  no  I/O) 

0 

'strio  ' 

= 

surf,  sain  trend  I/O  unit  (0  no  I/O) 

0 

'icvio  ' 

= 

ice  coverage  I/O  unit  (0  no  I/O) 

0 

'ithio  ' 

= 

ice  thickness  I/O  unit  (0  no  I/O) 

0 

'ictio  ' 

= 

ice  temperature  I/O  unit  (0  no  I/O) 

23 

'sshio  ' 

= 

sea  surf,  height  I/O  unit  (0  no  I/O) 
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0 

'  bsf io 

= 

baro .  strmfn. 

I/O 

unit 

(0 

no 

I/O) 

0 

' uvmio 

= 

mix.  lay. 

u-vel . 

I/O 

unit 

(0 

no 

I/O) 

0 

' vvmio 

= 

mix.  lay. 

v-vel . 

I/O 

unit 

(0 

no 

I/O) 

0 

' spmio 

= 

mix.  lay. 

speed 

I/O 

unit 

(0 

no 

I/O) 

0 

'  bltio 

= 

bnd.  lay. 

thick . 

I/O 

unit 

(0 

no 

I/O) 

0 

'mltio 

= 

mix.  lay. 

thick . 

I/O 

unit 

(0 

no 

I/O) 

0 

'  sstio 

= 

mix.  lay. 

temp . 

I/O 

unit 

(0 

no 

I/O) 

0 

' sssio 

= 

mix.  lay. 

sain . 

I/O 

unit 

(0 

no 

I/O) 

0 

' ssdio 

= 

mix.  lay. 

dens . 

I/O 

unit 

(0 

no 

I/O) 

1 

'kf 

= 

first  output  layer  (  = 

=  0  end 

output;  < 

layer 

1 

0 

#) 

'kl 

' uvlio 

= 

last  output  layer 
layer  k  u-vel .  I/O 

unit 

(0 

no 

I/O) 

0 

' vvlio 

= 

layer  k 

v-vel . 

I/O 

unit 

(0 

no 

I/O) 

0 

'  splio 

= 

layer  k 

speed . 

I/O 

unit 

(0 

no 

I/O) 

0 

'  inf io 

= 

layer  k 

i . dep . 

I/O 

unit 

(0 

no 

I/O) 

0 

'  thkio 

= 

layer  k 

thick . 

I/O 

unit 

(0 

no 

I/O) 

0 

'  temio 

= 

layer  k 

temp 

I/O 

unit 

(0 

no 

I/O) 

0 

'  salio 

= 

layer  k 

sain . 

I/O 

unit 

(0 

no 

I/O) 

0 

'  tthio 

= 

layer  k 

dens , 

I/O 

unit 

(0 

no 

I/O) 

0 

'  sfnio 

= 

layer  k 

strmfn . 

I/O 

unit 

(0 

no 

I/O) 

0 

'kf 

= 

first  output  layer  (  = 

=  0  end 

output;  < 

layer 

EOF 


#) 


} 

# 

#  -  extract  3D  fields  from  a  s 

# 

ncdfSd  0  { 

# 

#  -  interpolate  to  3D  z-levels 

#  -  and  output  to  netCDF.  Uses 

# 

CDF033=$ {NCDIR} /u3d"_" $ { DOMAIN } " 
CDF034=$ {NCDIR} /v3d"_"${ DOMAIN } " 
CDF035=$ {NCDIR} /t3d"_" $ { DOMAIN } " 
CDF036=$ {NCDIR} /s3d"_" $ { DOMAIN} " 
# 

touch  $CDF033  $CDF034  $CDF035 
/bin/rm  $CDF033  $CDF034  $CDF035 
# 


gle  HYCOM  archive  file 


from  a  single  HYCOM  archive  file 
linear  interpolation. 

" $ { E } "_"$ { dtg } $ { tau } . nc 
" $ { E } $ { dtg } $ { tau } . nc 
" $ { E } $ { dtg } $ { tau } . nc 
" $ { E } $ { dtg } $ { tau } . nc 

$CDF036 

$CDF036 


export  CDF  TITLE="HYCOM  $DOMAIN  $EXPT" _ 

[export  CDF  INST="Naval  Research  Laboratory"! 

export  CDF033=$CDF033 
export  CDF034=$CDF034 
export  CDF035=$CDF035 
export  CDF036=$CDF036 
$AWarc/archv2ncdf 3z  <<  EOF 
$ { HYCOMINPUTFILE } 

MERSEA 

000  'iexpt  '  =  experiment  number  xlO  (000=from  archive  file) 

3  'yrflag'  =  days  in  year  flag  ( 0=360 J1 6 , 1=3 66 J1 6, 2=366 JOl , 3- 


actual ) 

$NX  'idm 

$NY  ' j  dm 

32  'kdm 

34.0  ' thbase ' 


longitudinal  array  size 
latitudinal  array  size 
number  of  layers 
reference  density  (sigma  units) 
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0  'smooth'  =  smooth  the  layered  fields  (0=F,1=T) 

1  'iorign'  =  i-origin  of  plotted  subregion 

1  'jorign'  =  j-origin  of  plotted  subregion 

0  'idmp  '  =  i-extent  of  plotted  subregion  (<=idm;  0 

idm) 

0  'jdmp  '  =  j-extent  of  plotted  subregion  (<=jdm;  0 

j  dm) 

1  ' itype  '  =  interpolation  type  (0=sample, l=linear) 


50  'kz 

» 

=  number  of  levels 

0.000 

*  z 

f 

= 

sample 

depth 

0.500 

*  z 

» 

= 

sample 

depth 

1.083 

*  z 

f 

= 

sample 

depth 

1.762 

*  z 

f 

= 

sample 

depth 

2.554 

*  z 

f 

= 

sample 

depth 

3 .477 

*  z 

f 

= 

sample 

depth 

4.552 

*  z 

f 

= 

sample 

depth 

5.806 

*  z 

f 

= 

sample 

depth 

7.268 

*  z 

» 

= 

sample 

depth 

8 . 971 

*  z 

f 

= 

sample 

depth 

10 . 957 

*  z 

» 

= 

sample 

depth 

13.271 

*  z 

f 

= 

sample 

depth 

15.968 

*  z 

» 

= 

sample 

depth 

19.112 

*  z 

f 

= 

sample 

depth 

22 .777 

*  z 

f 

= 

sample 

depth 

27 .049 

*  z 

f 

= 

sample 

depth 

32 . 027 

*  z 

» 

= 

sample 

depth 

37 . 831 

*  z 

f 

= 

sample 

depth 

44.595 

*  z 

f 

= 

sample 

depth 

52 .479 

*  z 

f 

= 

sample 

depth 

61.669 

*  z 

f 

= 

sample 

depth 

72.380 

*  z 

» 

= 

sample 

depth 

84 .865 

*  z 

f 

= 

sample 

depth 

99.418 

*  z 

» 

= 

sample 

depth 

116.380 

*  z 

f 

= 

sample 

depth 

136.151 

*  z 

» 

= 

sample 

depth 

159.195 

*  z 

f 

= 

sample 

depth 

186.055 

*  z 

f 

= 

sample 

depth 

217.363 

*  z 

f 

= 

sample 

depth 

253 . 855 

*  z 

f 

= 

sample 

depth 

296.390 

*  z 

f 

= 

sample 

depth 

345.968 

*  z 

» 

= 

sample 

depth 

403 .755 

*  z 

» 

= 

sample 

depth 

471 .110 

*  z 

f 

= 

sample 

depth 

549.619 

*  z 

» 

= 

sample 

depth 

641 .128 

*  z 

f 

= 

sample 

depth 

747 .789 

*  z 

» 

= 

sample 

depth 

872 .111 

*  z 

f 

= 

sample 

depth 

1017 .019 

*  z 

» 

= 

sample 

depth 

1185 . 922 

*  z 

f 

= 

sample 

depth 

1382 .793 

*  z 

» 

= 

sample 

depth 

1612.263 

*  z 

f 

= 

sample 

depth 

1879.729 

*  z 

» 

= 

sample 

depth 

2191 .483 

*  z 

f 

= 

sample 

depth 

2554 .860 

*  z 

f 

= 

sample 

depth 

2978.406 

*  z 

» 

= 

sample 

depth 

3472 .000 

*  z 

f 

= 

sample 

depth 

4047 .000 

*  z 

» 

= 

sample 

depth 

4718.000 

*  z 

f 

= 

sample 

depth 

implies 

implies 
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5500.000  'z  '  =  sample  depth 


0 

'botio  ' 

= 

bathymetry 

I/O  unit 

(0 

no 

I/O) 

0 

'mltio  ' 

= 

mix . 1 . thk . 

I/O  unit 

(0 

no 

I/O) 

0 

' tempml ' 

= 

temperature 

jump  across 

mixed-layer  (degC,  0  no 

I/O) 

0 

' densml ' 

= 

density 

jump  across 

mixed-layer  (kg/m3,  0  no 

I/O) 

0 

'  inf io ' 

= 

interface  depths  I/O 

unit 

(0  no  I/O) 

0 

'wvlio  ' 

= 

w-velocity 

I/O  unit 

(0 

no 

I/O) 

33 

'uvlio  ' 

= 

u-velocity 

I/O  unit 

(0 

no 

I/O) 

34 

'vvlio  ' 

= 

v-velocity 

I/O  unit 

(0 

no 

I/O) 

0 

'splio  ' 

= 

speed 

I/O  unit 

(0 

no 

I/O) 

35 

'temio  ' 

= 

temperature 

I/O  unit 

(0 

no 

I/O) 

36 

'salio  ' 

= 

salinity 

I/O  unit 

(0 

no 

I/O) 

0 

'tthio  ' 

= 

density 

I/O  unit 

(0 

no 

I/O) 

EOF 

} 

################################################## 

E=$Eenv 

GLOBAL=$AREA 

DOMAINDIRHOME=$ARCH_area 

EXPTDIRSCR=$WRK_dir 

DOMAIN=hycom 

TOPOVER=$Tenv 

TOPO=$Topo_area 

# 

#  -  hycom  archive  output  which  are  input  files  to  this  program 

# 

OUTPUT=$ {NC_area } 
mkdir  -p  $OUTPUT 
INPUT=$ { DOMAINDIRHOME } 

# 

#  -  output  directory  for  output  (netcdf) 

# 

NCDIR=$OUTPUT 

mkdir  -p  $NCDIR  #  should  already  be  there 
LOGS=$OUTPUT/logs 
mkdir  -p  $LOGS 
# 

cd  $EXPTDIRSCR 

# 

touch  regional . depth . a  regional . depth . b 
/bin/rm  -f  regional . depth . [ab] 

In  -sf  $TOPO/depth  ${GLOBAL}  ${TOPOVER}.a  regional . depth . a 
In  -sf  $TOPO/depth  ${GLOBAL}  ${TOPOVER}.b  regional . depth . b 
# 

touch  regional . grid. a  regional . grid. b 
/bin/rm  -f  regional . grid. [ab] 

In  -sf  $TOPO/regional . grid . a  regional . grid . a 
In  -sf  $TOPO/regional . grid . b  regional . grid . b 
# 

alias  jul2cal=$MODASbin/doy2idtgl 

# 

#  -  loop  thru  files  (1  pero  snapshot  in  julian  days) 

# 

for  run  in  ncdf2d  ncdfSd;  do 
yyyy=$year 
day=$ jull 
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dend=$ jul2 

while  [  $day  -le  $dend  ]  ;  do 


fday='echo  $day  |  awk  ' {printf  "%3.3d",  $1}'' 
dtg='jul2cal  $day  $yyyy' 
echo  $dtg 

#  The  outer  nests  have  the  experiment  in  front  of  the  file 
hr="00  06  12  18" 
hr="00" 

for  tau  in  $hr 
do 

HYCOMINPUTFILE=${ INPUT} /archv_$ {GLOBAL} . $ { yyyy }_$ { f day }_$ { tau } . a 

if  [  !  -e  $HYCOMINPUTFILE  ] ;  then 

echo  "file  not  found  $HYCOMINPUTFILE"  >$LOGS/$ {yyyy}  ${fday}.err 
#  keep  processing 
f  i 

$run  $HYCOMINPUTFILE 

done 

day='echo  $day+l  |  be' 
done 
done 


3.0  RELO  NCOM  TEST  CASES  FOR  THE  NEW  JERSEY  COAST 

This  section  briefly  describes  the  sensitivity  of  Relo  NCOM  to  the  vertical  resolution  of  the 
boundary  conditions  obtained  from  GOFS  V3.0.  A  regional  Relo  NCOM  is  configured  for  the 
area  off  the  New  Jersey  coast  bounded  by  the  region  76°-68°W,  37.5°-43.5°N.  The  horizontal 
resolution  is  3  km  in  both  longitude  and  latitude  and  it  is  configured  with  40,  50  and  100  vertical 
levels  (40z,  50z  and  lOOz,  respectively).  In  each  configuration,  it  receives  GOFS  V3.0  BCs  with 
the  same  vertical  resolution  and  the  40z  case  uses  the  same  pre-defined  levels  noted  in  section 
2.0.  It  uses  wind  forcing  from  the  18  km  resolution  Coupled  Ocean  Atmosphere  Mesoscale 
Prediction  System  (COAMPS)  and  thermal  forcing  from  the  0.5°  Navy  Operational  Global 
Atmospheric  Prediction  System  (NOGAPS).  It  includes  river  discharge  (Barron  and  Smedstad, 
2002)  and  barotropic  tidal  forcing  from  Egbert  and  Erofeeva  (2002)  is  specified  at  the  open 
boundaries.  The  hindcasts  span  -'2  months  (22  May  2009  -  31  July  2009)  and  observations  are 
assimilated  via  NCODA  with  the  exception  of  profile  data.  These  are  withheld  to  be  used  as  an 
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independent  validation  dataset.  GOFS  V3.0  BCs  are  provided  to  the  inner  model  every  6  hours 
along  the  outermost  grid  row/column.  NCODA  is  turned  off  within  a  10  point  buffer  zone  along 
each  sidewall.  The  analysis  period  spans  the  entire  month  of  June  2009.  Figure  1  shows  the 
location  of  the  independent  observations. 


43°N 

42°N 

41°N 

40°N 

39°N 

38°N 


130.0 
27.5 
25.0 
E  22.5 

E  20.0 
E  17.5 
E  15.0 


76°W  75°W  74°W  73°W  72°W  71  °W  70°W  69°W  68°W 


Figure  1:  Sea  surface  temperature  (°C)  for  15  June  2009  from  New  Jersey  Relo  NCOM.  The 
triangles  indicate  the  location  of  the  independent  observations.  All  are  Argo  profiles  with  the 
exception  of  one  moored  buoy  in  the  northeast  corner  of  the  domain  that  provided  observations 
most  days  of  the  analysis  month. 


Similar  to  the  GOFS  V3.0  Phase  II  VTR,  a  temperature  versus  depth  error  analysis  was 
performed  using  independent  profiles  and  is  shown  in  Figure  2.  In  the  upper  200  m  of  the  water 
column,  50z  Relo  NCOM  has  reduced  mean  error  and  root  mean  square  error  (RMSE)  than  the 
40z  and  100  z  configurations. 
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Mean  Error  (°C)  RMSE  (°C)  Meon  Temp  (°C) 


Figure  2:  Temperature  (°C)  versus  depth  error  analysis  in  the  upper  200  mfor  New  Jersey  Relo 
NCOM  for  June  2009  using  53  unassimilated  profiles  depicted  in  Figure  1.  The  left  panel  is 
mean  error  (bias),  the  middle  panel  is  RMSE  and  the  right  panel  an  average  temperature  of  all 
profiles.  The  gray  curve  is  the  observations  and  Relo  NCOM  configured  with  40  levels  (blue),  50 
levels  (red)  and  100  levels  (black). 

An  error  analysis  of  the  acoustieal  proxy  variables  is  also  eomputed  and  shown  in  Table 
1 .  The  observational  profiles  did  not  all  contain  enough  information  to  compute  each  variable, 
and  so  the  number  of  independent  profiles  in  also  listed  in  column  2.  Note,  that  the  analysis  for 
SLD  is  only  based  on  eight  profiles.  The  mean  of  each  variable  is  listed  in  column  3.  The  region 
is  defined  by  shallow  mixed  and  sonic  layer  depths  during  this  time  frame.  Again,  50z  Relo 
NCOM  is  the  configuration  that  produces  the  lowest  overall  error. 

The  difference  in  these  results  is  a  consequence  of  the  interpolation  methodology 
employed  here.  The  GOFS  V3.0  native  grid  files  are  linearly  interpolated  onto  the  vertical  levels 
of  the  netCDF  files,  and  these  are  further  interpolated  onto  the  sigma-z  levels  of  the  Relo  NCOM 
domain  using  the  Piecewise  Cubic  Hermite  Interpolating  Polynomial  (PCHIP).  Linear 
interpolation  provides  line  segments  with  sharp  comers  at  the  data  points,  while  higher  order 
interpolations  are  made  for  matching  not  only  the  data  values  but  also  the  slopes  and  concavities 
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of  each  interpolating  segment.  The  PCHIP  method  modifies  the  derivative  values  in  the  Hermite 
representation  in  order  to  eliminate  the  bumps  and  wiggles  that  are  frequently  corrupting  the 
fields  interpolated  by  cubic  splines  or  Akiwa  algorithms  (Fritsch  and  Carlson,  1980).  Therefore, 
the  number  of  z  levels  in  the  netCDF  files  must  be  chosen  to  preserve  the  accuracy  of  GOFS 
V3.0  hybrid  vertical  coordinates,  yet  without  over  constraining  the  interpolated  variables  to  a 
high  number  of  linear  segments  in  such  a  way  that  the  PCHIP  algorithm  may  perform  in  an  area 
of  strong  gradients  (i.e.  the  thermocline).  In  our  example  the  40  levels  are  not  sufficient  to 
provide  a  correct  parameterization  of  the  GOFS  V3.0  thermocline  and  the  100  levels  are  over 
constraining  the  derivative  values  for  an  efficient  application  of  the  PCHIP  method. 


Table  1:  Acoustical  proxy  error  analysis  for  New  Jersey  Relo  NCOM  configured  for  40,  50  and 
100  levels  for  June  2009  using  the  number  of  unassimilated  profiles  noted  in  column  2.  Those 
cells  highlighted  in  green  have  the  lowest  bias  or  error  of  the  three  configurations. 


Variable 

#  of 
profiles 

Ohs. 

Mean 

40z 

50z 

lOOz 

Mean 

Error 

RMSE 

Mean 

Error 

RMSE 

Mean 

Error 

RMSE 

MLD 

52 

4.9  m 

0.0  m 

5.6  m 

-0.4  m 

5.1  m 

-0.6  m 

5.2  m 

SLD 

8 

30.0  m 

6.1  m 

9.7  m 

-1.9m 

2.4  m 

-1.9  m 

2.4  m 

BLG 

49 

5.3 

m/s/100  ft 

0.7 

m/s/100  ft 

2.3 

m/s/ 100  ft 

0.6 

m/s/100  ft 

2.1 

m/s/100  ft 

0.8 

m/s/ 100  ft 

2.2 

m/s/ 100  ft 

DSC 

34 

278.9  m 

-18.9  m 

74.8  m 

-17.8  m 

95.1  m 

-26.4  m 

88.6  m 

4.0  SUMMARY 

This  report  documents  the  two  methodologies  for  extracting  GOFS  V3.0  boundary 
conditions  for  use  in  Relo  NCOM.  In  the  first  method  netCDF  files  from  the  global  system  are 
on  a  constant  .08°  latitude/longitude  grid  and  have  been  vertically  remapped  to  40  pre-defmed 
levels.  The  second  methodology  allows  more  flexibility  it  uses  the  GOFS  V3.0  archive  files  on 
the  native  hybrid  vertical  grid  that  can  be  remapped  in  the  vertical  to  a  user  defined  set  of  z- 
levels. 
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A  version  of  Relo  NCOM  was  configured  for  the  area  off  the  New  Jersey  coast  and 
configured  for  40,  50  and  100  levels,  receiving  outer  BCs  from  GOFS  V3.0  at  the  same  vertical 
resolution.  Using  independent  profile  observations,  error  analyses  of  temperature  as  a  function  of 
depth  and  acoustical  proxy  variables  indicate  the  configuration  with  50  levels  had  overall  slightly 
lower  error  than  the  other  configurations. 
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6.0  APPENDIX 

The  following  describes  an  example  of  how  to  process  GOFS  V3.0  archive  files  where 
output  exists  only  over  a  pre-defined  subregion  of  the  global  domain.  These  are  on  the  hybrid 
vertical  grid  and  can  be  created  as  the  system  runs  in  real-time.  They  are  designated  as  the 
*archt*  files  and  the  subregion  is  relative  to  the  i/j  grid  as  defined  in  the  archt. input  file.  A 
sequence  of  scrip ts/software  must  be  run  to  process  and  make  them  ready  for  the  next  step  of 
extracting  BCs  for  Relo  NCOM,  i.e.  section  2.2. 

If  the  subregion  is  new,  fdes  that  define  the  latitude-longitude  grid  and  topographic 
depths  must  be  created.  This  need  only  be  done  once.  As  mentioned  above,  the  subregion  is 
defined  in  archt.input  and  may  look  something  like: 

esloginl  92>  cat  archt.input 
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0 

1 

0386  1112  0314  0331 

where  the  last  line  defines  the  subregion  and  is  listed  as:  first  i-eoordinate,  last  i-eoordinate,  first 
j-eoordinate,  and  last  j-eoordinate.  Again,  these  are  relative  to  the  global  domain.  If  the  EXACT 
latitude  and  longitude  of  the  lower  left  and  upper  right  eorners  of  the  subregion  are  known,  then 
prepO.make_ij.com  ean  be  used  to  extraet  the  i,j  eoordinates. 

davinci  98%  csh  prepO.make  ij.com 

USAGE  :  prepO.make  ij.com  -lonl  xx*  -lon2  xx*  -latl  xx* 

lat2  XX* 

EXAMPLE:  prepO.make  ij.com  -lonl  283.25  -lon2  292.75  -latl  35.60 
lat2_42.85 

davinci  99%  csh  prepO.make  ij.com  -lonl  105.0  -lon2  130.0  -latl  -30.0 
-lat2_-5 . 0 

i  386  j  1112  i dm  314  j dm  331 


Use  prepl.make_grid.com  to  create  the  subregion  file  that  defines  the  latitudes  and 


longitudes  and  prep2.make_depth.com  to  create  the  subregion  file  that  defines  the  topographic 


depths.  Note:  RRRrO.08  should  be  replaced  with  a  subregion  identifier  name  and  environmental 


variable  $T  is  the  topography  version  number  used  in  GOFS  V3.0. 


prep  1  .make_grid.com 

# ! /bin/csh 
# 

#  Purpose:  Create  the  Sub-Regional  &  Sub-Depth  files. 

#  Input  Files:  The  Global  Regional  &  Depth  files 

# 

#  Executable  for  Area  :  sub  grid 

#  Executable  for  Depth:  isub  topog 

# 

#  Version:  2.0  2011 

#  By  :  D.S.  Franklin,  QinetiQ-NA 

# 

################################################################### 

# 

#  Gather  these  variables  to  the  necessary  points 

# 

set  debug  =  0 

# 

foreach  arg  (  ${argv}  ) 
switch  (  ${arg}  ) 
case  -lonl  *  : 

setenv  LONl  'echo  ${arg}  |  sed  ' s/-lonl_// ' ' 
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breaksw 
case  -lon2  *  : 

setenv  L0N2  'echo  ${arg}  |  sed  ' s/-lon2_// ' ' 
breaksw 
case  -latl_*  : 

setenv  LATl  'echo  ${arg}  |  sed  ' s/-latl_// ' ' 
breaksw 
case  -lat2_*  : 

setenv  LAT2  'echo  ${arg}  |  sed  ' s/-lat2_// ' ' 


breaksw 
case  -i  *  : 
setenv  11 

echo  ${arg}  | 

sed  ' 

s/-i_// ' ' 

breaksw 
case  -j  *  : 
setenv  J1 

echo  ${arg}  | 

sed  ' 

s/-j_// ' ' 

breaksw 
case  -idm  *  ; 
setenv  IDM 

'echo  ${arg} 

1  sed 

'  s/-idm  // ' ' 

breaksw 
case  -jdm  *  : 
setenv  JDM 

'echo  ${arg} 

1  sed 

'  s / - j  dm  / /  '  ' 

breaksw 
case  -d  : 

set  debug  = 
breaksw 
default : 

^  1 

echo  "USAGE  :  'basename  $0'  -latl  x*  -lonl  x*  -lat2  x* 

lon2  X*" 

echo  "USAGE  :  or" 

echo  "USAGE  :  'basename  $0'  -i  *  -j  *  -idm  *  -jdm  *" 
exit  1 
breaksw 
endsw 
end 

# 

#  Global  Region  is  currently  4500  x  3298 

# 


if  ( 

!  $?R  ) 

setenv  R 

"gLBaOVOl" 

#  Original  Global  Region 

if  ( 

!  $?S  ) 

setenv  S 

"RRRrO  .  OOl" 

#  New  SubRegional  Name 

if  ( 

!  $?T  ) 

setenv  T 

"|09|" 

#  Bathymetry 

if  ( 

!  $?SRC 

)  setenv 

SRC 

1” 

'wallcraf /hvcom/ALL/bin" 

if  ( 

!  $?L 

)  setenv 

L 

pwd' 

if  ( 

!  $?T1 

)  setenv 

T1 

"$ {HOME 

; } /hycom/ $ {R} / topo" 

if  ( 

!  $?T2 

)  setenv 

T2 

" $ { HOME } /hycom/ $ { S } / topo " 

##################################################################### 

# 

#  No  need  to  change  anything  from  below  this  point 

# 

##################################################################### 

# 

#  Make  sure  coordinates  exist 

# 

if  (  !  $?I1  I  I  !  $?J1  I  I  !  $?IDM  I  I  !  $?JDM  )  then 
if  (  $?LAT1  &&  $?LAT2  &&  $?L0N1  &&  $?LON2  )  then 

setenv  IJINFO  '/bin/csh  ${L}/make  IJ.com  -latl  ${LAT1} 
lat2  ${LAT2}  -lonl  ${L0N1}  -lon2  ${L0N2}'  “  “ 
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setenv  II  'echo  ${IJINFO}  |  awk  '{  print  $2  }'' 

setenv  J1  'echo  ${IJINFO}  |  awk  '{  print  $4  }'' 

setenv  IDM  'echo  ${IJINFO}  |  awk  '{  print  $6 

setenv  JDM  'echo  ${IJINFO}  |  awk  '{  print  $8 

else 

echo  "USAGE  :  'basename  $0'  -latl  x*  -lonl  x*  -lat2  x*  -lon2  x*" 
echo  "USAGE  :  or" 

echo  "USAGE  :  'basename  $0'  -i  *  -j  *  -idm  *  -jdm  *" 
exit  2 
endif 
endif 

echo  "CREATING  DEPTH  FILE  for  ${S}  STARTING  AT  ${II}i 
COVERING  ${IDM}x${ JDM}  POINTS" 

# 

if  (  -d  ${TI}  )  then 
cd  ${TI} 
else 

echo  "NO  Directory  found:  ${T1}" 
exit  2 
endif 

if  (  !  -d  ${T2}  )  /bin/mkdir  -p  ${T2} 

# 

set  gridl  =  "${ T1 } /regional . grid" 
set  grid2  =  "${ T2 } /regional . grid" 

# 

#  Remove  any  previous  versions  of  the  subregional  files 

# 

if  (  -e  ${grid2}.a  )  /bin/rm  ${grid2}.a 
if  (  -e  ${grid2}.b  )  /bin/rm  ${grid2}.b 

# 

echo  "MAKING  GRID  FILE  for  Region  ${S}" 

${SRC}/. ./ subregion/ src/ sub_grid  <<E-o-D 
$ { gridl } . b 
$ { grid2 } . b 

${IDM}  'idm  '  =  longitudinal  array  size 

${JDM}  'jdm  '  =  latitudinal  array  size 

${11}  'irefi  '  =  longitudinal  input  reference  location 

$ { J1 }  'jrefi  '  =  latitudinal  input  reference  location 

1  ' iref o  '  =  longitudinal  output  reference  location 

1  'jrefo  '  =  latitudinal  output  reference  location 

E-o-D 

if  (  -e  ${grid2}.a  )  then 
echo  "CREATED  ${grid2}.a" 
else 

echo  "UNABLE  TO  CREATE  ${grid2}.a" 
exit  3 
endif 

# 

#  End  Program 

prep2.make_depth.com 

# ! /bin/csh 

#  SubRegDepth.com 

# 

#  Purpose:  Create  the  Sub-Regional  &  Sub-Depth  files. 

#  Input  Files:  The  Global  Regional  &  Depth  files 

# 

#  Executable  for  Area  :  sub  grid 
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#  Executable  for  Depth:  isub  topog 

# 

#  Version:  2.0  2011 

#  By  :  D.S.  Franklin,  QinetiQ-NA 

# 

################################################################### 

# 

#  Gather  these  variables  to  the  necessary  points 

# 


set  debug  =  0 


foreach  arg  (  ${argv}  ) 
switch  (  ${arg}  ) 
case  -lonl  *  : 


setenv  LONl  'echo  ${arg} 

1  sed 

'  s/-lonl  // ' ' 

breaksw 

case  -lon2  *  : 

setenv  LON2  'echo  ${arg} 

1  sed 

'  s/-lon2  // ' ' 

breaksw 

case  -latl  *  : 

setenv  LATl  'echo  ${arg} 

1  sed 

'  s/-latl  // ' ' 

breaksw 

case  -lat2  *  : 

setenv  LAT2  'echo  ${arg} 

1  sed 

' s/-lat2  // ' ' 

breaksw 

case  -i  *  : 

setenv  11  'echo  ${arg}  | 

sed  ' 

s/-i_// ' ' 

breaksw 

case  -j  *  : 

setenv  J1  'echo  ${arg}  | 

sed  ' 

s/-j_// ' ' 

breaksw 

case  -idm  *  : 

setenv  IDM  'echo  ${arg} 

1  sed 

' s/-idm  // ' ' 

breaksw 

case  -jdm  *  : 

setenv  JDM  'echo  ${arg} 

1  sed 

' s / - j  dm  / /  '  ' 

breaksw 

case  -d  : 

set  debug  =  1 

breaksw 

default : 

echo  "USAGE  :  'basename  $0' 

11 

-latl  X*  -1 

2^ 

echo  "USAGE  :  or" 

echo  "USAGE  :  'basename 

$0'  -i 

*  -j  *  -idm  * 

exit  1 

breaksw 

endsw 

end 


# 

#  Global  Region  is  currently  4500  x  3298 

# 


# 


if  ( 
if  ( 
if  ( 

!  $?R  ) 

!  $?S  ) 

!  $?T  ) 

setenv  R 
setenv  S 
setenv  T 

"SLBaO  .081" 
"RRRrO  .00;" 
"0  9" 

#  Original  Global 

#  New  SubRegional 

#  Bathymetry 

Region 

Name 

if  ( 

!  $?SRC 

)  setenv 

SRC  "~wallcraf /hvcom/ALL/bin 

II 

if  ( 

!  $?L 

)  setenv 

L  'pwd ' 
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if  (  !  $?T1  )  setenv  T1  " $ { HOME } /hycom/$ {R} /topo" 
if  (  !  $?T2  )  setenv  T2  "$ {HOME } /hycom/$ { S } /topo" 

# 

##################################################################### 

# 

#  No  need  to  change  anything  from  below  this  point 

# 

##################################################################### 

# 

#  Make  sure  coordinates  exist 

# 

if  (  !  $?I1  I  I  !  $?J1  I  I  !  $?IDM  I  I  !  $?JDM  )  then 
if  (  $?LAT1  &&  $?LAT2  &&  $?L0N1  &&  $?L0N2  )  then 

setenv  IJINFO  Vbin/csh  ${L}/make  IJ.com  -latl  ${LAT1} 
lat2_${LAT2}  -lonl_$ { LONl }  -lon2_$ { L0N2 } '  “  “ 

setenv  II  'echo  ${IJINFO}  |  awk  '{  print  $2  }'' 

setenv  J1  'echo  ${ IJINFO}  |  awk  '{  print  $4  }'' 

setenv  IDM  'echo  ${IJINFO}  |  awk  '{  print  $6  }'' 

setenv  JDM  'echo  ${IJINFO}  |  awk  '{  print  $8  }'' 

else 

echo  "USAGE  :  'basename  $0'  -latl  x*  -lonl  x*  -lat2  x*  -lon2  x*" 
echo  "USAGE  :  or" 

echo  "USAGE  :  'basename  $0'  -i  *  -j  *  -idm  *  -jdm  *" 
exit  2 
endif 
endif 

echo  "CREATING  DEPTH  FILE  for  ${S}  STARTING  AT  ${II}i  ${Jl}j 
COVERING  ${IDM}x${ JDM}  POINTS" 

# 

if  (  -d  ${TI}  )  then 
cd  ${TI} 
else 

echo  "NO  Directory  found:  ${T1}" 
exit  2 
endif 

if  (  !  -d  ${T2}  )  /bin/mkdir  -p  ${T2} 

# 

set  depthl  =  " $ { T1 } /depth_$ { R}_$ { T } " 
set  depth2  =  " $ { T2 } /depth_$ { S }_$ { T } " 

# 

#  Remove  any  previous  versions  of  the  subregional  files 

# 

if  (  -e  ${depth2}.a  )  /bin/rm  $ { depth2 } . a 
if  (  -e  ${depth2}.b  )  /bin/rm  ${depth2}.b 

# 

#  Create  the  depth  file.  This  requires  a  special  depth  file  to  take 
into 

#  account  the  tiled  regions.  The  resulting  global  file  will  be  large 
as  all 

#  the  cells  are  written  and  filled. 

# 

echo  "MAKING  DEPTH  FILE  for  ${S}" 

${SRC}/. . /subregion/src/isub_topog  <<E-o-D 
$ { depthl } . b 
$ { depth2 } . b 

${ depthl}  subregioned  to  ${S}  via  isub_topog 
${IDM}  'idm  '  =  longitudinal  array  size 

${JDM}  'jdm  '  =  latitudinal  array  size 
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${11}  'irefi  '  =  longitudinal  input  reference  location 

$ { J1 }  'jrefi  '  =  latitudinal  input  reference  location 

1  ' iref o  '  =  longitudinal  output  reference  location 

1  'jrefo  '  =  latitudinal  output  reference  location 

1  'ijgrd  '  =  integer  scale  factor  between  input  and  output 

grids 
E-o-D 

if  (  -e  ${depth2}.a  )  then 
echo  "CREATED  ${depth2}.a" 
else 

echo  "UNABLE  TO  CREATE  ${depth2}.a" 
exit  4 
endif 

# 

#  End  Program 

A  global  HYCOM  bathymetry  file  specifically  for  the  tiled  region  is  necessary  to  create 
the  subregion  HYCOM  archive  files.  This  file  should  exist  as 
/u/home/${user}/hycom/GLBa0.08/topo/depth_GLBa0.08_tileRRR.[ab].  However  this  needs  a 
combined  global  HYCOM  archive  file  of  the  tiled  regions  that  can  be  created  using  the  script 
prep3.make_tile.com.  The  IDM  and  JDM  are  for  the  entire  global  domain.  This  resultant 
bathymetry  file  puts  land  everywhere  the  HYCOM  file  has  a  p-grid  data  void.  If  this  is  not  done, 
there  will  be  the  following  error; 


error  -  wrong  bathymetry  for  this  archive 
number  of  topo  sea  mismatches  =  <num> 
number  of  topo  land  mismatches  =  <num> 

prep3.make_tile.com 

# ! /bin/csh 

# 

#  Program:  make  tile.com 

#  Create  a  masked  tiled  Global  Depth  file 

#  -  bathyetry  for  expt  74.2  subregion  archive. 

#  Since  directory  structure  changes  for  the  archt  files, 

#  pass  directory  information  via  command  line 

# 

#  Version:  2.0  2011 

#  By  :  D.S.  Franklin,  QinetiQ-NA 

#  Modified  Script  from  Alan  Wallcraft  (NRL) 

# 

################################################ 

# 

setenv  R  "GLBaO.08" 

setenv  V  "RRRrO.08" 

setenv  U  "echo  ${V}  |  cut  -cl-3' 

setenv  E  "|185|" 
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setenv  X  'echo  ${E}  |  awk  '{  printf ( " %04 . If \n" ,  $1  *  0.1  ) 

setenv  T  "[0~^" 
setenv  Y03^  |1 0 8| 

setenv  YY  'echo  ${Y03}  |  awk  '{  print  $1  +  1900 

setenv  A  "h" 
setenv  0  "${YY}  244" 
setenv  G  "01" 

# 


setenv  SRC  |"~wallcraf /hycom/ALL/bin^ 
setenv  T1  " $ { HOME } /hYcom/$ { R} /topo" 
setenv  T2  "${ HOME } /hycom/ $ {V} /topo" 
setenv  TMPF  "tmp.lis" 
setenv  OD  " $ { T1 } /depth_$ { R}_$ { T } " 
setenv  GD  "depth_$ {R}_$ {T }_tile$ {U} " 

setenv  S  "/scr/${user} /hycom/ $ {R} / expt_$ {X } /data/tart_$ { Y03 } $ { A} " 

setenv  L  "${ HOME } /hycom/ $ {V} /subregion" 
setenv  GF  " $ { S } /$ { E }_archt . $ { 0 }_$ { G } " 

# 

cd  ${S} 

# 

#  Check  for  existance  of  a  a  full  days  worth  of  global  archive  data 

#  Create  if  missing. 

# 

if  (  !  -e  ${GF}.a  )  then 

# 

#  Link  topo  files  for  global  processing 

# 

if  (  -e  . /regional . grid. a  )  /bin/rm  . /regional . grid. a 
if  (  -e  . /regional . grid. b  )  /bin/rm  . /regional . grid. b 
/bin/cp  $ { HOME } /hycom/$ {R} /topo/regional . grid . a  . 

/bin/cp  ${ HOME } /hycom/$ {R} /topo/regional . grid . b  . 

Create  a  a  full  days  worth  of  global  archive  data 


# 

# 

# 

awk  -f  $ { L } /subreg . awk  o=${0}  g=${G}  r=${R}  e=${E}  x=${X}  t=${T} 
s="${S}"  \ 

$ { L } / stepl . make  global.com  >! 

$ { E }_make_global_$ {0 }_$ { G} .com 

csh  $ { E }_make_global_$ {0 }_$ { G} . com  >&! 

$ { E }_make_global_$ {0 }_$ { G} .log 
endif 

# 

#  -  bathyetry  for  expt  ${X}  subregion  archive. 

# 

if  (  -e  ${GF}.a  )  then 

if  (  -e  ${GD}.a  )  /bin/mv  ${GD}.a++ 
if  (  -e  ${GD}.b  )  /bin/mv  ${GD}.b++ 

setenv  IDM  ' grep  idm  ${ T1 } /regional . grid . b  |  awk  '{  print  $1 
setenv  JDM  'grep  jdm  ${ T1 } /regional . grid . b  |  awk  '{  print  $1 

# 

if  (  -e  ${TMPF}  )  /bin/rm  -f  ${TMPF} 

${SRC}/hycom_mask  ${0D}.a  ${GF}.a  ${IDM}  ${JDM}  ${GD}.a  >  ${TMPF} 
if  (  !  ${status}  )  then 

echo  "Failed  to  create  ${GD}.a,  exiting" 
exit  5 
endif 
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#  -  uses  min, max  from  hycom  mask 

# 

set  min  =  ' grep  min  ${TMPF}  |  awk  ' {  print  $4  } ' ' 
set  max  =  ' grep  max  ${TMPF}  |  awk  ' {  print  $5  } ' ' 
cat  «E-o-D  >!  ${GD}.b 
from  depth_$ {R}_$ { T } 

MASKED  EVERYWHERE  EXCEPT  FOR  EXPT  ${X}  SUBREGION  ARCHIVE  LOCATIONS 
based  on  smith  and  sandwell  Imin 

Depths  greater  than  6500.0  scaled  by  0.20 
E-o-D 

printf  "min, max  depth  =  %10.3f  %9.3f\n"  $min  $max  »  ${GD}.b 

# 

cat  $ { CD } . b 

/bin/mv  -f  ${GD}.[ab]  ${T1}/ 

else 

echo  "ERROR  :  Missing  Tiled  Global  Archive  File  ${GF}" 
endif 

# 

#  End  Program 

In  the  main  subregional  topo  directory  a  new  file  of  the  type  "RegionIJ.csh"  is  required 
to  have  the  IDM,  JDM,  IREF,  JREF  and  file  size  information  in  the  C-Shell  environmental  setup 
for  future  sourcing.  This  is  created  by  running  prep4.make_RegionIJ.com. 


prep4.make_RegionIJ.com 

# ! /bin/csh 

setenv  R  "feLBaO . 08 " 

setenv  V  "^RRrO.08" 

setenv  W  p -wallcraf /hycom/ALL/bin/ "| 
setenv  I  "Region  ${V}  IJ.csh" 

# 

#  Determine  if  the  file  already  exists, 

#  if  so  back  it  up  into  a  unique  name  before  recreating. 

# 

if  (  -e  ${ HOME } /hycom/ $ {V} /topo/$ { I }  )  then 

/bin/mv  -f  $ { HOME } /hycom/$ { V} /topo/$ { I } 

$ { HOME } /hycom/ $ {V} /topo/ $ { I } . ' date  +"%Y%m%d%H%M" ' 
endif 
# 

#  From  the  Subregion  regional  grid  file,  find  the  IJ  domain 

# 

setenv  IDM  'grep  idm  $ {HOME } /hycom/$ {V} /topo/regional . grid . b  |  awk  '{ 
print  $1  } ' ' 

setenv  JDM  'grep  jdm  ${ HOME } /hycom/$ {V} /topo/regional . grid .  b  |  awk  '{ 
print  $1  } ' ' 

# 

#  Use  the  Subregion  regional  grid  file  to  determine 

#  the  start  of  the  region  in  global  lonlat  coordinates . 

#  Use  this  lonlat  to  determine  the  reference  point  on  the  global 
regional  grid. 

# 
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set  11  =  '${W}/hycom  ij21onlat 

$ { HOME } /hycom/ $ {V} /topo/ regional . grid,  a ' 
set  11  =  'echo  ${11}  |  sed  -e  's/E//'  -e  's/N//'  -e  '  s/\ (  . *\) W\ ( . *\) /- 

\1\2/'  -e  's/\(.*\)  \(.*\)S\(.* 

\) /\1  -\2/  '  ' 


set  11  =  ( ' $ {W} /hycom_lonlat2i j 

$ { HOME } /hycom/ $ { R} /topo/ regional . grid. a ' ) 


# 

#  Create  the  regional  IJ  file 

# 

echo  '  setenv  I DM  "'${IDM} 

echo  '  setenv  JDM  "'${JDM} 

echo  '  setenv  IREFI  "'$11[1] 

echo  '  setenv  JREFI  "'$11  [2] 


T  II  I 
I  II  I 
I  II  I 
I  II  I 


>!  $ { HOME } /hycom/$ { V} /topo/ $ { I } 
>>  $ { HOME } /hycom/$ { V} /topo/ $ { I } 
>>  $ { HOME } /hycom/$ { V} /topo/ $ { I } 
>>  $ { HOME } /hycom/$ { V} /topo/ $ { I } 


${11} 


Creating  the  tiled  HYCOM  archive  file  uses  the  steps  in  the  script 
stepO.submitwrap.com.  Modify  this  with  the  experiment  number,  the  subregions,  etc.  This 
queue  submission  wrapper  script  creates  a  global  HYCOM  archive  file  per  hour 
(stepl.make_global.com),  each  regional  subfde  per  hour  (step2.isubreghrly.com),  and  finally 
will  tar  up  the  resultant  subregional  hourly  files  into  daily  archives  (step3.tarhrfiles.com). 

All  scripts  are  created  using  the  HYCOM  technique  of  an  awk  script,  "subreg.awk",  to 
modify  the  resultant  scripts  as  needed.  The  user  should  only  have  to  modify 
stepO.submitwrap.com,  and  no  changes  should  be  required  in  stepl.make_global.com, 
step2.isubreghrly.com,  and  step3.tarhrfiles.com. 


step0.submitwrap  .com 

# ! /bin/csh 

#PBS  -N  XXX 

#PBS  -j  oe 

#PBS  -o  XXX. log 

#PBS  -1  walltime=6 : 00 : 00 

#PBS  -W  umask=027 

#PBS  -A  ACCTNUM 

#PBS  -q  standard 

# 

#  Wrapper  script  to  run  all  the  processing 

# 

#  Version:  2.0  2011 

#  By  :  D.S.  Franklin,  QinetiQ-NA 

#  Based  on  Scripts  by  E.J. Metzger,  NRL 
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# 

set  echo 
set  verbose 
# 

setenv  OS  'uname' 
switch  ($0S) 
case  ' SunOS ' : 
case  ' Linux ' : 
which  yod 

if  (!  $status)  then 

setenv  OS  XT3 
endif 

which  aprun 

if  (!  $status)  then 

C  -  XT4  or  XT5 

setenv  OS  XT4 
set  APRUN= ' aprun  -n  1' 
set  APRUNPLT= ' aprun  -n  1  ' 
endif 
breaksw 
case  'AIX' : 

set  APRUN=' ' 
set  APRUNPLT= ' ' 
breaksw 
default : 
set  APRUN=' ' 
set  APRUNPLT=' ' 


endsw 


# 

#  -  convert  tiled  to  standard  archives 

# 

#  Set  the  following  switches 

#  --  1  processing 

#  --  0  skip 

#  --  global  -  create  hourly  global  archive  file 

#  --  isub  -  create  hourly  subregional  file 

#  --  tarNrm  -  tar  up  subregional  files  and  remove  global  and 

subregional  files 

# 

set  global  =  1 
set  isub  =  1 
set  tarNrm  =  0 


# 

setenv  R  GLBaO.08 
setenv  E 
setenv  X 

setenv  T  _ 

setenv  VS  "RRRr 0.081" 
setenv  I  "Region  ${VS} 


IJ . csh" 


# 

#  -  environment  variables  defining  the  year  and  part 

# 

setenv  Y03  |108| 

setenv  YY  'echo  ${Y03}  |  awk  '{  print  $1  +  1900 

setenv  A  h 

setenv  GLB  HOURS  "|00  01  02  03  04  05  06  07  08  09  10  11  12  13  14  15  1 
117  18  19  20  21  22  23l" 
setenv  DTGl  |214| 
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42 


setenv  DTG2  |2  4  4| 
setenv  Z  "pTol^" 

# 

#  -  environment  variables  defining  the  working  directories 

# 

setenv  S  /scr/$ {user } /hycom/$ {R} /expt_$ {X } /data/tart_$ { Y03 } $ { A} 
setenv  L  $ { HOME } /hyc2ncom/subregion 
# 

date  +"STATUS  :  %D  %T  :  Using  arguments 

date  +"STATUS  :  %D  %T  :  REGIONS  =  ${VS}" 

date  +" STATUS  :  %D  %T  :  GLB_HOURS  =  $ { GLB_HOURS } " 

if  (  !  $ {global}  )  date  +"WARNING:  %D  %T  :  SKIPPING  GLOBAL  Processing" 
if  (  !  ${isub}  )  date  +"WARNING:  %D  %T  ;  SKIPPING  REGIONAL 

Processing  unless  missing" 

# 

while  (  'echo  ${DTG1}  |  sed  's/^0*//''  <=  'echo  ${DTG2}  |  sed 

's/^0*//''  ) 


# 

# 

# 


# 

# 

# 


if  (  'echo  ${DTG1}  |  wc  -c'  >  8  )  then 

If  Ordinal  day  is  from  DTGl :  Output  is  in  form  YEAR  ORD 

setenv  Y  'echo  ${DTG1}  |  cut  -cl-4' 

setenv  MTH  'echo  ${DTG1}  |  cut  -c5-6' 

setenv  DAY  'echo  ${DTG1}  |  cut  -c7-8' 

setenv  0  ' $ { L } /ymd2doy . csh  ${Y}  ${MTH}  ${DAY}' 

else 

setenv  0  $ {YY}_$ {DTGl } 
endif 

date  +"STATUS  :  %D  %T  :  Processing  for  ${0}" 

GLOBAL  PROCESSING 
if  (  $ {global}  )  then 


cd  ${S} 

# 

#  Link  topo  files  for  global  processing 

# 

if  (  -e  . /regional . grid. a  )  /bin/rm  . /regional . grid. a 
if  (  -e  . /regional . grid. b  )  /bin/rm  . /regional . grid. b 
/bin/cp  $ { HOME } /hycom/$ {R} /topo/regional . grid . a  . 

/bin/cp  ${ HOME } /hycom/$ {R} /topo/regional . grid . b  . 

# 

#  Create  a  full  days  worth  of  global  data 

# 

foreach  G  (  ${GLB_HOURS}  ) 

awk  -f  $ { L } /subreg . awk  o=${0}  g=${G}  r=${R}  e=${E}  x=${X}  t=${T} 
z=${Z}  s="${S}"  \ 

$ { L } /stepl . make  global.com  >! 

$ { E }_make_global_$ {0 }_$ { G} .com 

csh  $ { E }_make_global_$ {0 }_$ { G} . com  >&! 

$ { E }_make_global_$ {0 }_$ { G} . log  & 

# 

#  Set  automatically  the  size  of  the  file  for  the  required  region 

# 

if  (  -e  $ { E }_archt . $ {0 }_$ { G} . a  )  then 

set  sz  =  'Is  -1  $ {E }_archt . $ {0 }_$ { G} . a  |  awk  '{  print  $5  }'' 
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=4=  =4=  =«=  =«=  =4=  =«=  =«=  <r>  <r>  rt  -c/>  =t^=  =«=  =«= 


if  (  "${sz}"  >  "${Z}"  )  setenv  Z  ${sz} 
endif 
end 
wait 

Sometimes  the  global  file  is  not  created,  rerun  if  missing 

foreach  G  (  ${GLB_HOURS}  ) 

if  (  !  -e  $ {E }_archt . $ {0 }_$ { G} . a  )  then 

date  +"WARNING:  %D  %T  :  First  pass  did  not  create 
{E }_archt . $ {0 }_$ { G} . a,  trying  again" 

awk  -f  $  { L }  /  sutjreg .  awk  o=${0}  g=${G}  r=${R}  e=${E}  x=${X} 

=${T}  z=${Z}  s="${S}"  \ 

$ { L } /stepl .make_global . com  >!  $ { E }_make_global_$ {0 }_$ { G} . com 

csh  $ { E }_make_global_$ {0 }_$ { G} . com  >&! 

{ E }_make_global_$ {0 }_$ { G} .log 

if  (  !  -e  $ {E }_archt . $ {0 }_$ { G} . a  )  then 
date  +"ERR0R  :  %D  %T  :  Unable  to  create 
{E}_archt.${0}_${G} .a" 
endif 
endif 
end 

echo  "FINISHED  Creating  Global  Data  file" 
endif 

REGIONAL  PROCESSING 

if  (  ${isub}  )  then 
foreach  V  (  ${VS}  ) 

setenv  U  'echo  ${V}  |  cut  -cI-3' 

cd  ${S} 

if  (  -e  $ { HOME } /hycom/ $ { V} /topo/ $ { I }  )  then 

Source  specific  dimensions  of  region 

/bin/cp  ${ HOME } /hycom/ $ {V} /topo/ ${ I }  . 
source  $ { HOME } /hycom/ $ { V} / topo/ $ { I } 

Create  regional  archive  file 

foreach  G  (  ${GLB_H0URS}  ) 

Check  that  a  tile  file  has  been  created 


set  tfile  =  $ { HOME } /hycom/$ { R} /topo/depth_$ { R}_$ { T }_tile$ { U} 
if  (  !  -e  ${tfile}.a  )  then 

awk  -f  $ { L } /subreg . awk  o=${0}  g=${G}  r=${R}  e=${E}  x=${X} 

t=${T}  \ 

v=${V}  idm=${IDM}  jdm=${JDM}  ir=${IREFI} 
jr=${JREFI}  s="${S}"  \ 

$ { L } /prep3 . make  tile.com  >!  ${E}  make  tile  ${V}.com 
/bin/csh  ${E}  make  tile  ${V}.com  >&! 

$ { E }_make_tile_$ {V} .log 
endif 

/bin/cp  -p  ${tfile}.  [ab]  . 

if  (  -e  $ { E }_archt . $ {0 }_$ { G} . a  )  then 

awk  -f  ${ L } /subreg . awk  o=${0}  g=${G}  r=${R}  e=${E}  x=${X} 

t=${T}  \ 
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v=${V}  idm=${IDM}  jdm=${JDM}  ir=${IREFI} 
jr=${JREFI}  s="${S}"  \ 

$ { L } / step2 . isubreghrly . com  >! 

$ { E }_isubreghrly_$ {0}-${G}_${V} . com 

/bin/csh  $ {E }_isubreghrly_$ { 0 } -$ { G }_$ { V} . com  >&! 

$ { E }_isubreghrly_$ {0}-${G}_${V}.log  & 
else 

date  +"ERROR  :  %D  %T  :  Global  file  $ {E }_archt . $ {0 }_$ { G} . a 
missing,  can  not  create  subregion" 
endif 
end 
wait 
else 

date  +"ERR0R  :  %D  %T  :  Missing  file 
$ { HOME } /hycom/ $ { V} / topo/ $ { I } " 

date  +"ERR0R  :  %D  %T  :  Creating  the  file  necessary  to  have 
the  IDMxJDM  and  length  of  region,  Rerun" 

awk  -f  $ { L } /subreg . awk  r=${R}  v=${V} 

$ {L} /prep4 .make  RegionIJ.com  >!  ${E}  make  RegionIJ.com 

/bin/csh  ${E}  make  RegionIJ.com  >& !  ${E}  make  RegionIJ.log 
endif 
end 
endif 

# 

#  At  end  of  the  day,  after  the  last  hourly  file: 

# 

#  Remove  hourly  files  if  directed 

# 

if  (  ${tarNrm}  )  then 
foreach  V  (  ${VS}  ) 

awk  -f  ${ L } /subreg . awk  o=${0}  r=${R}  e=${E}  x=${X}  t=${T}  \ 

v=${V}  gh="${GLB_HOURS}"  s="${S}"  \ 

$ { L } /step3 . tarhrf iles . com  >! 

$ { E }_tarhrf iles_$ { 0 }_$ { V} . com 

/bin/csh  $ {E }_tarhrf iles_$ {0 }_$ { V} . com  >&! 

$ { E }_tarhrf iles_$ { 0 }_$ { V} . log  & 
end 
endif 

# 

date  +"STATUS  :  %D  %T  :  End  of  Processing  for  ${0}" 

date  +" STATUS  :  %D  %T  :  =====================================" 

# 

#  Update  the  day  by  your  favorite  program 

# 

if  (  'echo  ${DTG1}  |  sed  's/^0*//''  >  19000000  )  then 

perl  $ { L } /dtgadj . pi  ${DTG1}  +dl 
else 

setenv  DTGl  'expr  ${DTG1}  +  1' 
endif 

date  +"STATUS  :  %D  %T  :  New  date  is  ${DTG1}" 
end 
wait 


step  1  .make_global.com 

# ! /bin/csh 

# 

#  Version:  2.0  2011 

#  By  :  D.S.  Franklin,  QinetiQ-NA 
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#  Based  on  Scripts  by  E.J. Metzger,  NRL 

# 

set  echo 
set  verbose 
set  gstatus  =  0 
set  gsize  =  0 
# 

#  -  convert  tiled  to  standard  global  archive  file 

# 

if  (  !  $?APRUN  )  setenv  APRUN  "" 

setenv  R  GLBaO.08 
setenv  E  742 

setenv  X  'echo  ${E}  |  awk  '{  printf ( "%04 . If \n" ,  $1*0.1  )  }' 

setenv  T  09 
setenv  0  2007  001 
setenv  G  00 
# 

#  -  environment  variables  defining  the  working  directories 

#  -  tart_2007083118_2007090100 

# 

setenv  S  " /scr/$ {user } /hycom/$ {R} /expt_$ {X } /data" 

setenv  Z  "10000000000" _ 

setenv  SRC  |"  ~wallcraf /hycom/ALL/bin^ 

# 


date 

+ 

"STATUS 

%D 

%T 

Using 

arguments 

date 

+ 

"STATUS 

%D 

%T 

YEAR  ORD 

=  ${0} 

date 

+ 

"STATUS 

%D 

%T 

HOUR 

=  ${G} 

date 

+ 

"STATUS 

%D 

%T 

Min 

File  Size 

=  ${Z} 

# 

######################################################### 

# 

#  Create  GLOBAL  file  from  the  tiled  files 

# 

######################################################### 

# 

if  (  !  -d  ${S}  )  mkdir  -p  ${S} 

# 

#  Run  over  the  number  of  hours  in  a  day 

# 

cd  ${S} 

date  +"STATUS  :  %D  %T  :  Processing  Day  ${0}  for  Hour  ${G}" 

if  (  !  -e  archt2archv  )  /bin/cp  $ { SRC } / . . /archive/src/archt2archv  . 

# 

#  Link  topo  files  for  global  processing 

# 

if  (  !  -e  regional . grid. a  )  /bin/cp 
$ {HOME } /hycom/ $ {R} /topo/ regional . grid. a 
if  (  !  -e  regional . grid. b  )  /bin/cp 
$ { HOME } /hycom/ $ { R} /topo/ regional . grid. b 
# 

set  gfile  =  " $ { E }_archt . $ {0 }_$ { G} " 

# 

#  Remove  previous  versions  of  the  file 

# 

/bin/rm  -f  ./${ gfile }. [ab]  >& !  /dev/null 
echo  " . / $ { gf ile } . a"  >!  ${gfile}.tmp 

# 

date  +"STATUS  :  %D  %T  :  Searching  for  ./???? ?/archt .${ 0 }_${ G }. A" 
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/bin/ls  -1  . / ? ?? ??/archt . $ {0 }_$ { G} . A  »  ${gfile}.tmp 
set  files  =  'wc  -1  ${gfile}.tmp  |  awk  '{  print  $1 


date  +"STATUS  :  %D  %T  :  Running  archt2archv  with  $ {files}  files 
${0}  ${G}  with  $ { gf ile } . tmp" 

${APRUN}  . /archt2archv  <  ${gfile}.tmp 
set  gstatus  =  ${status} 
set  gsize  =  0 

if  (  -e  ${gfile}.a  )  set  gsize  =  'Is  -1  ./${gfile}.a  |  awk  '{ 
print  $5  } ' ' 

if  (  ${gsize}  <  ${Z}  )  then 

if  (  -e  ${gfile}.a  )  /bin/rm  -f  $ { gf ile } . [ab] 
set  gstatus  =  5 
endif 

if  (  $ {gstatus}  !=  0  )  date  +"ERR0R  :  %D  %T  :  Problem  creating 
global  files  for  ${0}  ${G}" 


for 


step2.isubreghrly.com 

# ! /bin/csh 

# 

#  Version:  2.0  2011 

#  By  :  D.S.  Franklin,  QinetiQ-NA 

#  Based  on  Scripts  by  E.J. Metzger,  NRL 

# 

set  echo 
set  verbose 
# 

#  -  convert  tiled  to  standard  archive  file 

# 

if  (  !  $?APRUN  )  setenv  APRUN  "" 

setenv  R  GLBaO.08 
setenv  E  742 

setenv  X  'echo  ${E}  |  awk  '{  printf ( "%04 . If \n" ,  $1*0.1  )  }' 

setenv  T  09 

setenv  0  2007_001 

setenv  G  00 

# 

#  Set  regional  variables 

# 

setenv  V  "ITFtO.08" 

setenv  U  'echo  ${V}  |  cut  -cl-3' 

setenv  IBM  "XXX" 

setenv  JDM  "XXX" 

setenv  IR  "XXX" 

setenv  JR  "XXXX" 


# 

#  -  environment  variables  defining  the  working  directories 

# 

setenv  S  "/scr/$ {user} /hvcom/$ {R} /expt  $ {X} /data" 
setenv  W  |"~wallcraf /hycoin/ALL/bin"^ 

# 

setenv  T1  " $ { HOME } /hycom/$ { R} /topo" 
setenv  T2  "${ HOME } /hycom/$ {V} /topo" 

# 


date  +"STATUS  : 
date  +" STATUS  : 
date  +" STATUS  : 

# 

#  Create  regional 
for  it 

# 

date  +" STATUS  : 


%D  %T  :  Using  arguments 
%D  %T  :  YEAR_ORD  =  ${0}" 

%D  %T  :  REGION  =  ${V}" 

hourly  files,  this  needs  a  specific  tiled  depth 


%D  %T  :  Creating  regional  files  for  ${0}  ${G}" 


file 


created 
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# 

#  Change  to  regional  directory,  link  needed  files,  if  necessary 

# 

if  (  !  -d  ${S}/${V}  )  mkdir  -p  ${S}/${V} 

cd  ${S}/${V} 

if  (  !  -e  isubregion  )  /bin/cp  $ { W} / . . /subregion/src/isubregion  . 

if  (  !  -e  regional . grid . a  )  /bin/cp  ${ T1 } /regional . grid . a  . 

if  (  !  -e  regional . grid . b  )  /bin/cp  ${ T1 } /regional . grid . b  . 

# 

set  GD  =  "depth_$ { R}_$ { T }_tile$ { U } " 

if  (  !  -e  ${GD}.a  )  /bin/cp  $ { T1 } /$ { GD } . a  . 

if  (  !  -e  ${GD}.b  )  /bin/cp  $ {T1 } /$ {GD} .b  . 

# 

set  RD  =  "depth_$ { V}_$ { T } " 

if  (  !  -e  ${RD}.a  )  /bin/cp  $ { T2 } /$ { RD } . a  . 

if  (  !  -e  ${RD}.b  )  /bin/cp  $ {T2 } /$ {RD} .b  . 

# 

#  Set  input  and  output  name,  removing  any  previous  versions  of  the  output 

# 

set  input_name  =  " $ { S } /$ { E }_archt . $ { 0 }_$ { G} " 
set  output_name  =  "$ {E}_archt_$ {V} . $ {0}_$ {G} " 

/bin/rm  $ { output_name } . [ab]  >& !  /dev/null 

# 

#  Create  the  subregion  hourly  file 

# 

date  +"STATUS  :  %D  %T  :  Running  isubregion  for  ${V}  ${0}  ${G}" 

${APRUN}  ./isubregion  <<E-o-D 
$ { input_name } . b 
${GD}  .b 

$ { output_name } . b 
$ { RD } . b 

${R}  interpolated  to  ${V} 

${IDM}  ' idm  '  =  longitudinal  array  size 

${JDM}  'jdm  '  =  latitudinal  array  size 

${IR}  'irefi  '  =  longitudinal  input  reference  location 
${JR}  'jrefi  '  =  latitudinal  input  reference  location 
1  ' irefo  '  =  longitudinal  output  reference  location 

1  'jrefo  '  =  latitudinal  output  reference  location 

1  'ijgrd  '  =  integer  scale  factor  between  input  and  output  grids 

0  'iceflg'  =  ice  in  output  archive  flag  (0=none, l=energy  loan  model) 

0  'smooth'  =  smooth  interface  depths  (0=F,1=T) 

E-o-D 

step3.tarhrfiles.com 

# ! /bin/csh 

# 

#  At  end  of  the  day,  after  the  last  hourly  file: 

#  -  Collect  all  regions  into  tar  file 

#  -  Clean  up  hourly  files 

# 

#  Version:  2.0  2011 

#  By  :  D.S.  Franklin,  QinetiQ-NA 

#  Based  on  Scripts  by  E.J. Metzger,  NRL 

# 

set  echo 
set  verbose 
# 

setenv  R  GLBaO.08 
setenv  V  ITFtO.08 
setenv  E  742 
setenv  X  74.2 
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setenv  GH  "00  01  02  03  04  05  06  07  08  09  10  11  12  13  14  15  16  17  18  19 
20  21  22  23" 
setenv  HR  "" 
setenv  0  2007  001 
set  remove  =  0 
set  gstatus  =  0 
# 

#  -  environment  variables  defining  the  working  directories 

# 

setenv  S  / scr/ ${ user } /hycom/$ { R} /expt_$ { X} /data 
setenv  OUT_R  "${S}/${V}" 

# 

date  +"STATUS  :  %D  %T  :  Using  arguments 

date  +"STATUS  :  %D  %T  :  OUT_R  =  ${OUT_R}" 

date  +" STATUS  :  %D  %T  :  YEAR_ORD  =  ${0}" 

# 

#  Create  Tar  file  of  all  hours  in  a  day 

# 

cd  ${OUT_R} 

# 

set  tarfile  =  " $ {E }_archt_$ { V} . $ {0 } . tar" 
foreach  h  (  ${GH}  ) 

if  (  !  -e  " $ {E }_archt_$ { V} . $ {0 }_$ { h } . a"  )  then 
date  +"WARNING:  %D  %T  :  Missing  hourly  file 
${E}_archt_${V} . $ { 0 }_$ {h } . a" 
setenv  HR  "${HR}  ${h}" 
set  gstatus  =  115 
else 

if  (  -e  $ { S } /$ {E }_archt . $ { 0 }_$ {h } . a  )  then 

date  +"STATUS  :  %D  %T  :  Removing  Global  file 
${E}_archt.${0}_${h} .a" 

if  (  -e  $ { S } /$ {E }_archt . $ { 0 }_$ {h } . a  )  /bin/rm  -f 

${S}/${E}_archt.${0}_${h} .a 

if  (  -e  $ { S } /$ {E }_archt . $ { 0 }_$ {h } . b  )  /bin/rm  -f 

${S}/${E}_archt.${0}_${h} .b 

if  (  -e  $ { S } /$ {E }_archt . $ { 0 }_$ {h } . tmp  )  /bin/rm  -f 
$ { S } /$ {E }_archt . $ { 0 }_$ {h } . tmp 
endif 
endif 
end 

if  (  !  ${gstatus}  | |  ${gstatus}  ==  "115"  )  then 

date  +"STATUS  :  %D  %T  :  Creating  tar  file  ${tarfile}" 
tar  -cvf  ${tarfile}  $ { E }_archt_$ { V} . $ { 0 }_? ? . [ab] 
if  (  ${gstatus}  !=  "115"  )  set  gstatus  =  ${status} 

# 

#  Remove  hourly  files  after  creation  of  the  tar  file 

# 

date  +"STATUS  :  %D  %T  :  Completed  Gtar  for  tarfile  ${tarfile}  DTG 
${0}  with  status  =  ${gstatus}" 

if  (  ${gstatus}  ==  0  &&  -e  ${tarfile}  )  then 
Is  -1  ${tarfile} 
if  (  ${ remove}  )  then 

date  +"STATUS  :  %D  %T  :  Removing  hourly  Regional  Files  for  ${V} 

${0}  :" 

/bin/rm  -f  $ {OUT_R} /$ { E }_archt_$ { V} . $ { 0 }_? ? . [ab] 
endif 

else  if  (  ${gstatus}  ==  "115"  )  then 

date  +"WARNING:  %D  %T  :  Missing  hourly  files  :  ${HR}" 
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else 

date  +"WARNING:  %D  %T  :  Unable  to  create  tar  file  ${tarfile}" 
endif 
else 

date  +"WARNING:  %D  %T  :  Unable  to  create  tar  file  ${tarfile}" 
endif 
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8.0  TABLE  OF  ACRONYMS 


BCs 

Boundary  Conditions 

BLG 

Below  Layer  Gradient 

DSC 

Deep  Sound  Channel 

DSRC 

DoD  Supereomputing  Research  Center 

GOFS 

Global  Ocean  Forecast  System 

HYCOM 

HYbrid  Coordinate  Ocean  Model 

MLD 

Mixed  Layer  Depth 

MODAS 

Modular  Ocean  Data  Analysis  System 

NAVOCEANO 

Naval  Oceanographic  Office 

NCODA 

Navy  Coupled  Ocean  Data  Assimilation 

NCOM 

Navy  Coastal  Ocean  Model 

NLOM 

NRL  Layered  Ocean  Model 

NRL 

Naval  Research  Laboratory 

SLD 

Sonic  Layer  Depth 

T 

Temperature 

VTR 

Validation  Test  Report 
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